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ABSTRACT 


This report is a manual to aid in using the STEEP32 code. STEEP32 
is the EXEC VIII version of the STEEP code (STEEP is an acronym for Shock 
Two-dimensional Eulerian Elastic Plastic) . The major steps in a STEEP32 
run are illustrated in a sample problem. There is a detailed discussion of 
the internal organization of the code including a description of each 
subroutine. 
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SECTION 1 


INTRODUCTION 


STEEP32 is a modification of the STEEP code initially developed 
under NASA contract NAS8-20235. Under the current contract STEEP has 
been made compatible with the EXEC VIII system at Marshall Space Flight 
Center. The name STEEP is an acronym for Shock Two-dimensional 
Eulerian Elastic Plastic. The suffix 32 was added since the EXEC VIII 
version of STEEP is restricted to using less than 32K words of storage. 

The development of the basic Eulerian method used in STEEP is 
described in Reference 1 and in Appendix D of Reference 2 . Subsequent 
to the initial development, the formulation was generalized to include 
strain rate effects on material properties. A discussion of the strain 
rate formulation and an analysis of scaling relationships in hypervelocity 
impacts appears in Reference 2 . An application of STEEP to the study of 
debris formed by hypervelocity impacts on thin plates is described in 
Reference 3. Comparisons between STEEP solutions and dynamically- 
instrumented experiments confirming the accuracy of the code are made 
in Reference 4. 

The documentation of STEEP 32 concentrates on the data processing 
since a thorough description of the mathematical formulation already 
appears in References 1 and 2. However, formulas in these documents 
are frequently cited to correlate the mathematics with the code. 
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SECTION 2 

ORGANIZATION OF STEEP 3 2 


The simulation of a STEEP32 problem takes place in three stages. 
These are initial conditions, integration, and edit. The body of computer 
programs associated with each of these operations is referred to a s a 
processor. 

The initial conditions, integration, and edit processors communicate 
to one another through information carried on magnetic tape. Each block of 
such information is called a restart file. Thus the initial conditions 
processor writes a restart file which in turn acts as part of the input to the 
integration processor. Each execution of the integration processor writes 
at least one restart file which in turn serves as data to subsequent integra- 
tion or edit runs . 

In the following three sections the sequence of operations within 
each of the main processors of STEEP32 is discussed. 


2 . 1 INITIAL CONDITIONS 

The initial conditions processor first uses subroutine DATA to read 
the namelist INITAL. This namelist contains all the data for initial 
conditions other than that associated with the boundary string. The next 
activity of initial conditions is the grid generation. During this phase of 
the computation subroutine GRIDVT uses the information in NETBLK to 
generate the grid variables in common block RCDATA. 

Subroutine GENBND is next used to generate the string information. 
The data associated with the string is in the namelist BDSTRG. Using the 
string information and grid coordinates subroutine MOVER calculates the 
area and volume data. Subroutine REGION completes the calculations by 
computing the initial values for the cell variables in the common block 
GRDATA . 

Initial conditions terminates by making a pass over the drum to 
produce print and restart files. The process is controlled by subroutine 
OUTPUT. 


2 . 2 INTEGRATION 

The first activity of the integration processor is RESTART. Following 
restart the namelist DATAL2 is read to modify the values that were on tape. 
This is done by subroutine DATA 2 . 
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The actual computations during integration take place in three stages. 
The first stage is controlled bv subroutine COMP1. This stage of the 
computation is sometimes referred to as the "no flow" portion of the integra- 
tion since only estimated quantities are calculated. COMP1 makes one 
pass over the drum. During this stage of the calculations the string points 
are moved. 

The second stage of the integration computations is the boundary 
string calculations. During this process the updated locations of the 
string points are used to obtain new areas and volumes. As with initial 
conditions the subroutine which controls these calculations is MOVER. 

The final stage of the computations is controlled by subroutine 
COMP2. During this stage the drum file created by COMP1 is used to 
perform the "flow" calculations and produce the final values of the cell 
variables of the current cycle. 

The following two charts illustrate the staging of calculations 
within COMF1 and COMP2. Note that COMP1 requires that there be five 
column buffers but COMP2 needs only four. 


2 . 3 EDIT 

The edit processor produces SC4020 plots in two stages. The first 
stage draws the background information. This consists of a grid, x and y 
labels and heading information. The background is controlled by subroutine 
PAPER. 


The second stage of the edit processor draws the foreground 
information. This process is controlled by subroutine FUNCT. During the 
second stage one pass is made over the scratch drum file. 

The creation of SC4020 edits from STEEP32 is also a function of the 
MSFC SC4020 library. Since this library is still in a state of flux, it is 
recommended that anyone using this part of STEEP32 should remain 
sensitive to changes being made by the systems staff. 
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STAGING OF CALCULATIONS COMP 2 
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SECTION 3 


AN ELEMENT BY ELEMENT DESCRIPTION OF STEEP32 


To the EXEC VIII system each of the main processors of STEEP32 is 
a distinct program file made up of individual elements . In EXEC VIII there 
are four types of elements: 

(1) Procedure elements . 

(2) Source language elements. 

(3) Relocatable binary elements. 

(4) Absolute binary elements . 

Since the last two types of elements are produced from the first two, they 
will not be discussed further. 

Procedure and source language elements can be either assembler, 
COBOL, or FORTRAN. There can also be COLLECTOR source elements. 
In.STEEP32 only three types of elements occur: 

(1) FORTRAN procedure elements. 

(2) FORTRAN source elements . 

(3) COLLECTOR source elements. 

The STEEP32 elements are described in the following pages according to 
their type . 


3 . 1 FORTRAN PROCEDURE ELEMENTS 

The FORTRAN procedures or PROC's currently in STEEP32 are used 
for two purposes. The first use is to define the values for the PARAMETERS. 
These PARAMETERS serve mainly to allocate space for larger arrays used by 
STEEP32. By changing the values of the PARAMETERS in these PROC's and 
recompiling the FORTRAN source code, it is possible to tailor the amount of 
core required to the needs of the particular problem being simulated. 

The second use of the procedures is to define frequently appearing 
common blocks. Using procedures in this fashion eliminates the need to 
have a copy of a common block for each subroutine using it. In addition 
to the common blocks, these procedures may also contain INTEGER, REAL, 
or LOGICAL statements for variables in the common block. It is important to 
note that a common block which appears in a procedure may also still be 
referenced explicitly in some of the older subroutines. 
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On the following pages each procedure is documented. The descrip- 
tion of the procedures containing common blocks is brief since the common 
blocks themselves are documented in a separate section. 



PROCEDURE PARAM 


PARAM contains the parameters defining the dimensions of a 
majority of the arrays in STEEP32. 


Parameter 

Description 

NUMROW 

The maximum number of rows a 
problem ca n have . 

NUMCOL 

The maximum number of columns 
a problem can have. 

NUMCEL 

The maximum number of cells. 

NUMGSC 

The number of columns in core that 
are used for scratch purposes. 
NUMGSC is currently set to 5 and 
this value should never change. 

NUMCUT 

The number of cell variables. 

NUMRDT 

The number of row variables in 
RC DATA . 

NUMCDT 

The number of column variables in 
RC DATA . 

NUMRTP 

The number of row temporary words . 

NUMCTP 

The number of column temporary 
words . 

NUMSTG 

The maximum number of boundary 
strings . 

NUMSTR 

The maximum number of string points 

NUMCUT 

The maximum number of cut cells. 

NUMCTY 

The maximum number of cell types. 

NUMSTY 

The maximum number of segments 
per cell. 

NUMHED 

The number of header words. 

NUMCOE 

The number of coefficients for 
artificial viscosity. 
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Procedure PARAM continued . 


Parameter 

Description 

NUMMSK 

The number of masks in dictionary 
word . 

NUMPCN 

The maximum number of PC ON 
constants . 

NUMMAT 

The number of materials. This must 
always be 1 . 

NUMCSG 

The maximum number of cuts per 
string segment. 

numreg 

The maximum number of regions 
used to define initial conditions. 

NUMXYS 

The maximum number of X, Y groups 
which may be used to specify grid. 
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/ 

PROCEDURE PARAMB 

PARAMB contains the parameters used to define the dimension of 


those arrays used exclusively for the 
Parameter 
NUMSCG 
NUMWDH 

NUMTIT 

NUMDEP 

NUMSPC 

NUMCON 


SC4020 edits. 

Description 

The maximum number of SC4020 graphs 

The maximum number of words for 
labels. 

The maximum number of words for 
title . 

The maximum number of deposition 
bins. 

The maximum number of contour values 
used in table. 

Maximum number of contour values 
used in labeling the contours. 
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PROCEDURE BOUCOM 


BOUCOM contains the common block BOUNDY and several type 
statements. Thus it contains all the variables associated with the string and 
the string points. 


PROCEDURE CONCOM 

CONCOM contains the common block CONSTS and several type 
statements for variables in CONSTS. In addition, CONCOM defines an array 
ALLCOM which is equivalenced to CYCLE. This is to make it easy for sub- 
routine FILE to write the contents of CONSTS onto the restart tape. 


PROCEDURE CUTCOM 

CUTCOM contains the common block CUTBLK. This holds all the 
information about partially filled cells. 


PROCEDURE GDACOM 

GDACOM contains all the scratch area for the cell variables. In 
addition to a type statement GDACOM also defines a number of arrays which 
are in turn equivalenced to an appropriate location of GDATA. This makes it 
possible to reference any of the cell variables by use of a single subscript. 


PROCEDURE MATCOM 

MATCOM contains the common block MATBLK and a type statement 
for the logical variable VONMIS. 


PROCEDURE RCDCOM 

RCDCOM contains the common block RCDATA. RCDCOM also 
defines several arrays which it equivalences to appropriate parts of 
variables in RCDATA. This allows for single subscripts in those subroutines 
which must have acess to the grid information. 


PROCEDURE SENCOM 

SENCOM contains the common block SENSWT which used to hold the 
values of the first six sense switches. All references to SENCOM are 
believed to have been removed. 
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PROCEDURE SUBCOM 

SUBCOM contains common block SUBSBK. This is all the subscript 
Information necessary to reference any of the cell or grid variables. 


PROCEDURE SHHCOM 

SHHCOM contains common block SHCOM which defines the strain 
hardening information. SHHCOM also includes the type statement for the 
logical variable SF. 


3.2 


FORTRAN SOURCE ELEMENTS 


The majority of elements in the STEEP32 code are FORTRAN V source 
elements. Several of these elements such as BET and IBET were originally 
written in the SLEUTH assembly language for the EXEC II system. Those 
elements were changed to FORTRAN to avoid the possibility of encountering 
incompatibilities between the EXEC II and EXEC VTII assemblers. Other 
elements such as GEDATE and GETIME were originally written in SLEUTH to 
extract special words from the EXEC II resident. These are now dummy 
FORTRAN subroutines. Routines serving similar purposes should be 
obtainable from the EXEC VIII system library or staff. 

The following pages document the FORTRAN source elements. Under 
the heading "Element Description" the code is related to the mathematics 
described in several publications. Each element's role in the data 
processing is also described. 

Under "Abort Conditions" the references to subroutine ABORT are 
described. The first word of the description appears in all capital letters. 
ABORT prints this word preceded by the words "ABORT CALLED BY" . when 
it is called. STEEP32 calls subroutine ABORT when it detects an illogical 
situation occurring. 


The section entitled "Messages to Printer" describes what effect 
each subroutine has on the print file. 
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Element Name: 


ABORT 


Element Type: FORTRAN Subroutine 


Entry Points: ABORT 

Called By: ARECEL, BDNORM, BDSTRS, GENBND, GRIDVT , INTRST, 

PUTGET, SETMAT, SETSEG , SETSID, SETSUB, SIO, VOLCEL, 
BDMOVE , COLPTS, SELDT, STRESS, GETCON, PLASFL, 
ISOMET, FXVSFY, MASSES, PLOTTR, CONTUR 


External References: None 

Common Blocks: IOCNTL, SUBSBK 


Calling Sequence: ANAMS 

Processor: I.C. 


Variable 

Description 

ANAME 

1 word alphanumeric description of tb t e type of call. 


Element Description: ABORT terminates via an illegal computed GO TO. 

This forces a dump if the "PMD"' card is present. 


Abort Conditions: None 


Message to Printer: ABORT CALLED BY ANAME 
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Element Name: ARECEL _ Element Type: FORTRAN FUNCTION 

Subroutine 

Entry Points: ARECEL 

Called By: PRNTIJ, VOLCEL, SETFOR 

External References: GETC ON , SAVSUB, GETIJ, ABORT 
Common Blocks: CONSTS, SUBSBK, CUTBLK, RCDATA 
Calling Sequence: L, MODE 


Variable 

Description 

(1) L 

(2) MODE 

Legal Subscript to Cell Dictionary. 

= 1 , A Regular Cell is Assumed, 

f 

- = 2, A Stress Cell is Assumed. 


Element Description: ARECEL calculates the area of a material lying 

in either a regular or stress cell. 


Abort Conditions: ARECEL - If the value in the cell dictionary is 

not a positive number. 


Message to Printer: None 
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Element Name: ARVOL 


Element Type: FORTRAN Subroutine 


Entry Points: ARVOL 

Called By: BDNORM, BDSTRS 

External References: ATAN2 

Common Blocks: None 

Calling Sequence: (XT, YT, IT, XS/YS, IS, JS, ISIDE, XSIDE, AREA, VOL) 


Variable 

Description 

IT 

Number of string points and intersections of the 
string with the boundary. 

XT, YT 

Coordinates of the string. 

XS, YS 

Location of the sides of the cell. 

IS, JS 

Subscripts of the cell . 

ISIDE, XSIDE 

Initial and final sides which are being cut by the 
string . 

AREA , 
VOL j 

Output variables. 


Element Description: ARVOL uses line integrals to calculate the area and 

volume of the material bounded by the string points 
lying in a cell. 


Abort Conditions: None 


Messages to Printer: None 
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Element Name: BDINIT 

Entry Points: BDINIT 

Called By: MOVER 


Element Type: FORTRAN Subroutine 


External References: GETCON, PUTGET, PUTCON, INTPSH 


Common Blocks: CONTS, BOUNDY 

Calling Sequence: None 


Variable 

Description 




Element Description: BDINIT initializes the space associated with the 

boundary push down list. In the dictionary it 
replaces time N information with time N + 1 
information. 


Abort Condition: None 

Message to Printer: None 
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Element Name: BDMOVE 


Element Type: FORTRAN Subroutine 


Entry Points: BDMOVE 

Called By: COLPTS 

External References: EXTRAP, EXIST, SETSUB, ABORT, MASSIJ 

Common Blocks: CONSTS, SUBSBK, GRDATA, RC DATA , IOCNTL 

Calling Sequence: $,IM, JM , LIM, LJM , XM, YM 


Variable 

Description 

$ 

• Nonstandard return. 

IM 
JM J 

I,J coordinates of regular cell that string point is in. 


I,J coordinates of lattice cell that string point is in. 

U 

X,Y coordinates of string point. 


Element Description: BDMOVE moves any string point which is found to 

exist in a column. The double interpolation scheme 
used is described on page 27 of Reference 1 , how- 
ever, the mass terms have been deleted from the 
equations 3.35. The nonstandard return is taken if 
no cell surrounding the string point exists. 


Abort Conditions: BDMV.l - The denominator in equations 3.35 

is zero. 

BDMV.2 - A string point is outside the extremum 
of the problem. 


Messages to Printer: None 
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Element Name: 


BDNORM 


Element Type: FORTRAN Subroutine 


Entry Points: BDNORM 

Called By: BDNORM 

External References: ENCLOS, CHOOSE, SETSUB, ABORT, SETSEG, 

INTRST, ARVOL, STARVL 

Common Blocks: SUBSBK, RCDATA, BOUNDY, N TRACT 

Calling Sequence: (SI, SJ, X, Y) 


Variable 

Description 

SJ, SI 
X,Y 

' 

Cell subscripts of cell being worked with. 

Initial intersection point coordinates. 


Element Description: BDNORM collects the string points in a regular 

cell. It calls INTRST to determine the intersection 
of the string with the cell sides. It then calls 
ARVOL to determing the area and volume of the 
regular cell. 


Abort Conditions: BDNORM - If the point (X,Y) does not lie on one 

of the four cell-sides. 


Messages to Printer: None 
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Element Name: BDPOST Element Type: FORTRAN Subroutine 

Entry Points: BDPOST 

Called By: MOVER 

External References: GETCON, SETSUB, SETSID, FILLER, SETMAT 

Common Blocks: RCDATA, BOUNDY, CONSTS, SUBSBK, SCANBK, 

CUTBLK 

Calling Sequence: None 


Variable 

Description 




Element Description: Subroutine BDPOST cleans up the remaining items at 

the termination of processing the boundary. A 
scan is made over the areas and volumes of the cells 
to see whether the cut areas and cut volumes are indeed 
very close to the total areas and total volumes, in 
which case. the total areas and total volumes are sub- 
stituted for cut areas and cut volumes. This is done 
for both regular cells and stress cells. An additional 
computation prohibits cut cell areas and volumes 
from exceeding the magnitude of the actual cells . 
Following this sequence of computations a series of 
calls to SETSID is made in order to set horizontal 
and vertical lattice segment types and cell segment 
types in the remainder of the dictionary. Finally a 
call to subroutine FILLER and a call to subroutine 
SETMAT are made to complete the dictionary. 


Abort Conditions: None 


Messages to Printer: None 
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Element Name: 


BDSTRS 


Element Type: FORTRAN Subroutine 


Entry Points: BDSTRS 

Called By: BONDPR 

External References: ENCLOS, CHOOSE, SETSUB, ABORT, SETSEG, INTRST 

ARVOL, STARVE 

Common Blocks: SUBSBK, RCDATA, BOUNDY 

Calling Sequence: SI, SJ, X,Y 


Variable 

Description 

SI,SJ 

X,Y 

Cell subscripts of cell being worked with. 
Initial intersection point coordinates. 


Element Description: Same as BDNORM only for stress cells. 

Abort Conditions: BDSTRS - If the point (X,Y) does not lie on one of 

four sides of the lattice cell. 

Messages to Printer: None 
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Element Name: BET Element Type: FORTRAN FUNCTION 

Entry Points: BET 

Called By: ENC LOS, LINES, TBFT02, LNEDIT, CONTUR 

External References: None 

Common Blocks: None l 

Calling Sequence: A, X, B, M 


Variable 

Description 

A 

First test argument. 

X 

Second test argument. 

B 

Third test argument. 

: m 

Dummy argument 


Element Description: Routine checks if X is between A and B. 

See description of IBET. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: BONDPR 


Element Type: FORTRAN Subroutine 


Entry Points: BONDPR 

Called By: MOVER 

External References: FOLTES, PACUNP, BDNORM, BDSTRS 

Common Blocks: CONSTS, BOUNDY, IOCNTL 

Calling Sequence: None 


Variable 

Description 




Element Description: BONDPR is the driving routine for boundary 

processing. It calls BDNORM and BDSTRS to 
calculate the areas and volumes of regular and 
stress cells respectively. FOLTES is called to 
check for a string fold condition. 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: BOUNST Element Type: FORTRAN Subroutine 

Entry Points: BOUNST 

Called By: SETFOR 

External References: None 

Common Blocks: GRDATA, SEGBLK 

Calling Sequence: L, A, B, C 


Variable 

Description 




Element Description: BOUNST determines the stresses at a point of cut 

in a regular cell. Currently BOUNST returns the 
stresses indexed by L. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: CBN DRY Element Type: FORTRAN LOGICAL 

FUNCTION 

Entry Points: CBNDRY 

Called By: PHASEA 

External References: GETCON 

Common Blocks: CONSTS 

Calling Sequence: L, MODE 


Variable 

Description 

L 

MODE 

Pointer to dictionary word. 

= 1; Segment is a horizontal, regular cell segment. 
= 2; Segment is a regular cell vertical segment. 

= 3; Horizontal stress segment. 

= 4; Vertical stress segment 


Element Description: CBNDRY tests for the existence of the segment type 

indicated by MODE. CBNDRY= .TRUE . if a segment 
of the type desired exists. 


Abort Conditions: None 


Messages to Printer: None 
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Element Name: CHOOSE 


Element Type: FORTRAN INTEGER 

FUNCTION Subroutine 

Entry Points: CHOOSE 

Called By: BDNORM , BDSTRS 

External References: None 

Common Blocks: None 

Calling Sequence: X, Y, I, J, XS, YS 


Variable 

Description 

X,Y 

Point coordinates. 

IJ 

Cell subscripts . 

XS, YS 

Sides of the cell. 


Element Description: CHOOSE selects the cell sides nearest to the 

point (X, Y) . 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: CINOUT 

Entry Points: CINOUT 

Called By: MAIN 2 


Element Type: FORTRAN Subroutine 


External References: 
Common Blocks: 
Calling Sequence: 
Processor: 


RESTIN, RESTOU 

NONE 

NONE 

INTEG 


Variable 


Description 


Element Description: CINOUT is a driver routine for the restart processor. 

It was originally inserted to force a particular 
allocation in the EXEC II version of STEEP 32. 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: COLD Element Type: FORTRAN Subroutine 

Entry Points: COLD 

Called By: MAIN 1 , MAIN2 , MAIN 3 

External References: SETSUB 

Common Blocks: CONSTS, SHCOM., CUTBLK, DELSBK , SEGBLK , 

SUBSBK, C BLOCK, RCDATA 

Calling Sequence: None 

Processor: I.C., INTEG, EDIT 



Element Description: Subroutine COLD initializes several variables 

defining the grid and strain hardening points. 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: COLPTS Element Type: FORTRAN Subroutine 

Entry Points: COLPTS 

Called By: PHASEA 

External References: PACUNP, MAS SI J, ABORT, BDMOVE 

Common Blocks: IOCNTL, RC DATA , BOUNDY, SUBSBK 

Calling Sequence: None 

Processor: INTEG 


Variable 

Description 

ICOL 

IT 

POINTER to column under consideration 

= 1 if ICOL is left-most column in problem 
= 2 if ICOL is an interior column 
= 3 if ICOL is right-most column in problem 


Element Description: Subroutine COLPTS moves all the string points in 

a given column. If a point is found to lie within 
the given column and has not already been moved 
then COLPTS calls BDMOVE to move the point. 


Abort Conditions: COLPTS - If a regular string point was found to 

have no cells around it with material. 


Messages to Printer: None 
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Element Name: 
Entry Points: 
Called By: 

External Reference 

Common Blocks: 
Calling Sequence: 
Processor: 


COMPUT Element Type: FORTRAN Subroutine 

COMPUT 

MAIN 2 

s: SELDT, COMP I, MOVER, COMP2 , STPOST 

CONSTS, SHCOM, FRSTBK, MATBLK 

$ / L 
INTEG 


Variable 

Description 

$ 

Nonstandard return which is taken when there is an 


interrupt in the normal flow of the integration 


computations. 

L 

Currently takes on the value 1. 


Element Description: COMPUT is the driving routine for the math portion 

of the integration processor. It: is described with 
the integration processor. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: COMP1 


Element Type: FORTRAN Subroutine 


Entry Points: COMP1 

Called By: COMPUT 

External References: SOPEN, PNTSET, SRITER, SREAD, PHASEA, HYPOEL, 

SC LOSE 

Common Blocks: CONSTS 

Calling Sequence: None 

Processor: INTEG 



Element Description: COMP1 performs the first phase of the integration 

calculations. COMP1 makes one pass over the 
scratch drum file. 

Abort Conditions: None 

Messages to Printer:. None 
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Element Name: - COMP2 


Element Type: FORTRAN Subroutine 


Entry Points: 


COMP2 


Called By: 


COMPUT 


External References: 

Common Blocks: 
Calling Sequence: 
Processor: 


SOPEN, SREAD, STRFLO, MASFLO, STRESS, 
SRITER, SC LOSE 

CON STS 

None 

INTEG 


Variable 


Description 


Element Description: Subroutine COMP2 performs the second phase of 

the math portion of the integration processor. 
COMP2 uses as inputs the scratch file created by 
COMP1 and in turn outputs an updated set of cell 
variables or drum. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: 

CONTLB 

Element Type: FORTRAN 

Block Data Routine 

Entry Points: 

NONE 



Called By: 

External References: None 

Common Blocks: SCCNTL 

Calling Sequence: None 

Processor: EDIT 


Variable 


Description 


Element Description: CONTLB initializes constants which are used only 

by the SC4020 processor. Since it is a BLOCK 
DATA routine it has no executable statements. 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: CONTRL Element Type: FORTRAN 

Block Data Routine 

Entry Points: None 

Called By: 

External References: None 

Common Blocks: MATBLK, DELSBK, SCANBK, RCDATA, RCTEMP , 

STJBSBK, CONSTS, IOCNTL, TOLBLK 

Calling Sequence: None 

Processor: INTEG, EDIT 



Element Description: CONTRL is a BLOCK DATA program which initializes 

various items in the common blocks to their default 
values . 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: 

CONTUR 

Element Type: FORTRAN Subroutine 

Entry Points: 

CONTUR 


Called By: 

FUNCT 


External References: GETCON 

VARIBL, 

, ABORT, PLTBND , SOPEN, SREAD, SETSUB, 
BET 

Common Blocks: 

SUBSBR, 

RCTEMP 

CONSTS, SEGBLK, RCDATA, SCCNTL, 

Calling Sequence: 

Q 


Processor: 

EDIT 



Variable 

Description 

Q 

The subscript to be used in the variables defining 
a particular plot,- 


Element Description: Subroutine CONTUR plots equal value contours 

of selected cell quantities on the SC4020. 

CONTUR makes one complete pass over the 
drum. 

Abort Conditions: CONT.l If NCONT(Q) is not between 1 and NUMCON 

CONT.2 If CONT YP(Q) is not between 1 and 7 

Messages to Printer: None 
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Element Name: DATA Element Type: FORTRAN Subroutine 

Entry Points: DATA 

Called By:. . DUMSUB,- MAIN3 

External References: TITLE, SETSUB 

Common Blocks: NETBLK, REGBLK, MATBLK, CONSTS, SHCOM, 

IOCNTL, PCNBLK, SUBSBK, RCDATA, BOUNDY 

Calling Sequence: None 

Processor: I.C., EDIT 

Variable Description 


Element Description: DATA is used by the initial conditions and edit 

processors to read the namelist INITAL. DATA 
also initializes the region data to zero and reads 
in additional PCON constants if requested. For 
the edit processor DATA is used merely to define 
the I/O units. 

Abort Conditions: None 

Messages to Printer: DATA prints an abbreviated version to INITAL 

called INITLO. If additional PCON constants 
are read in, they too are listed. 
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Element Name: 


DIVU 


Element Type: FORTRAN Subroutine 


Entry Points: 
Called By: 

External Reference 
Common Blocks: 
Calling Sequence: 
Processor: 


DIVU 

HYPOEL, VISC OS 
s: None 

DIVUBK, SUBSBK, CONSTS, RCDATA 

IT, JT, U, V 

INTEG 


Variable 

Description 

IT 

= 1 if left-most column 


= 2 if interior column 


= 3 if right -most column 

JT 

= 1 if bottom row 


= 2 interior row 


= 3 top row 

U 

Velocity component in X direction. 

V 

Velocity component in Y direction. 


Element Description: Subroutine DIVU calculates the velocity gradients 

across a cell surface. The results are passed 
back to the calling routine via DIVUBK. The 
equations used in DIVU are 3.32 and 3.33 in 
Reference 1 . 

t . . .. • • 

. Abort Conditions: None 

.. Messages to Printer: None 
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Element Name: DMSUB1 Element Type: FORTRAN Subroutine 

Entry Points: DMSUB1 

Called By: MAIN1 

External References: REGION, OUTPUT 

Common Blocks: GRDATA , CBLOCK, NETBLK, REGBLK, CONSTS, 

PCNBLK 

Calling Sequence: None 

Processor: I.C. 



Element Description: DMSUB1 is used by the initial conditions processor 

to call subroutines REGION and OUTPUT. DMSUB1 
is then used in the MAP segment STARTR which 
isolates several subroutines and common blocks 
in the overlay structure . 

Abort Conditions: None 

Messages to Printer: None 


38 






Element Name: DOTPLT Element Type: FORTRAN Subroutine 

Entry Points: DOTPLT 

Called By: CONTUR 

External References: None 

Common Blocks: None 

Calling Sequence: XII, YI1, XI2, YI2, ID1, INCTR, KKK 

Processor: EDIT 



Element Description: DOTPLT is called by CONTUR to draw dotted lines 

for the pressure contours. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: DUMSUB Element Type: FORTRAN Subroutine 

Entry Points: DUMSUB 

Called By: MAIN 1 

External References: DATA 

Common Blocks: GRDATA , CBLOCK, NETBLK, REGBLK, CONSTS, 

PCNBLK 

Calling Sequence: None 

Processor: I.C. 



Element Description: DUMSUB is a dummy subroutine left from the 

EXECII version of STEEP. Its sole function is to 
call DATA . 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: 

ENCLOS Element Type: 

FORTRAN Logical 
Function 

Entry Point: 

ENCLOS 


Called By: 

BDNORM, BDSTRS 


External References: BET 


Common Blocks: 

None 


Calling Sequence: 

(CX, CY, TX, TY, TI) 


Processor: 

I.C., INTEG 



Variable 

Description 

CX 

CY 

TX 

TY 

TI 

X coordinate of cell center. 

Y coordinate of cell center. 

Array containing X coordinates of string points. 
Array containing Y coordinates of string points. 
The number of string points. 

1 

; 


Element Description: ENCLOS determines if the point (CX, CY) lies 

inside of the material. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: 


EXIST 


Element Type: FORTRAN Subroutine 


Entry Points: 
Called By: 


EXIST 

PRNTIJ, BDMOVE, HYPOEL, SETDEL, SETFOR, STRESS, 
STRFLO , MASFLO, PLTSTR, VARIBL, VECTOR 


External References: GETCON 

Common Blocks: . CONSTS 

Calling Sequence: (L, M) 


Processor: 


I.C., INTEG, EDIT 


Variable 

Description 

L 

Subscript of cell dictionary word. 

M 

= 1 

Regular cell center at time N. 


= 2 

Lattice point at time N . 


= 3 

Lattice point at time N . 


= 4 

Lattice point at time N-l . 


= 5 

Lattice point at time N and N-l. 


= 6 

Regular cell center at time N . 


= 7 

Regular cell center at time N + 1. 


Element Description: EXISTS calls subroutine GETCON to extract 

existence information for the points indicated 
by the value of M. 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: 

EXTRAP 

Element Type: FORTRAN Logical 
Function 

Entry Points: 

EXTRAP 


Called By: 

BDMOVE 


External References: IBET 


Common Blocks: 

CON STS, SUBSBK 


Calling Sequence: 

UA f JA, M) 


Processor: 

INTEG 



Variable 

Description 

IA 

I offset 

JA 

J offset 

M 

= 1 Regular cell 

= 2 Lattice cell 


Element Description: EXTRAP determines the existence of the point 

(I+IA, J+JA) or (LI+IA, LJ+ JA) depending on the 
value of M . 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: -FCUT -- - Element Type: FORTRAN Function 

Subroutine 

Entry Points: FCUT 

External References: 

Common Blocks: None 

Calling Sequence: A, B, C , D # Y, R 

Processor: INTEG 


Variable 

Description 

A 

0 


XX 

B 

a 


yy 

C 

°TT 

D 

a 


xy 

Y 

Yield 

■ 

R 

a.. • a.,/2 
lj \Y 


Element Description: FCUT is used in conjunction with the Perzyna model 

documented in Section 2.2 of Reference 2. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: 

FILE 

Element Type: FORTRAN Subroutine 

Entry Points: 

FILE 


Called By: 

RESTOU, RESTIN 

External References: 

SETSUB, SRITER, SC LOSE, SOPEN, GETIME, 
GEDATE, SREAD, PRNTIJ 

Common Blocks: 


CONSTS, MATBLK, SUBSBK, BOUNDY, PCNBLK, 
SHCOM, C BLOCK, GRDATA , RCDATA, CUTBLK, 
IOCNTL 

Calling Sequence: 


$, $, MODE, MATCH 

Processor: 


I.C., INTEG, EDIT 


Variable 

Description 

$ 

Nonstandard return. 

$ 

Nonstandard return. 

MODE 

= 1 Next restart file read from tape. 


= 2 Restart file is written. 


= 3 Skip to next file . 

MATCH 

FILE number to MATCH. 


Element Description: Subroutine FILE reads or writes a selected restart 

file. The first nonstandard return is taken if an 
end of tape was encountered on a read. The second 
nonstandard return is taken when an illegal or 
unexpected record type is encountered. 

Abort Conditions: None 


Messages to Printer: After each record other than a cell variable record 

is read, FILE prints a message indicating the record 
type. Once the machine is completely refreshed 
FILE prints a message indicating this and identifying 
the logical unit, file number, and cycle number which 
is used. 
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Element Name: 


FILLER 


Element Type: FORTRAN Subroutine 


Entry Points: FILLER 

Called By: BDPOST 

External References: GETCON, SETSUB, PUTCON 

Common Blocks: SUBSBK, CONSTS 

Calling Sequence: None 

Processor: I.C., INTEG 


Variable Description 



Element Description: FILLER completes entries in the cell dictionary 

after the boundary has been processed. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: FOLTES Element Type: FORTRAN Logical 

Function 

Entry Points: FOLTES 

Called By: BONDPR 

External References: LINES 

Common Blocks: None 

Calling Sequence: N, X, Y 



Element Description: FOLTES tests for a fold among inputted segments of 

strings. Each segment is tested against non-adjacent 
segments through subroutine LINES. If an intersection 
is found, LINES takes the alternate return and FOLTES 
is set to true . 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: FORDIC 

Entry Points: FORDIC 

Called By: MOVER 

External References: SETSUB, PUTCON 

Common Blocks: CONSTS, SUBSBK 

Calling Sequence: None 

Processor: I.C.,INTEG 


Element Type: FORTRAN Subroutine 


Varia ble 


Description 


Element Description: 
Abort Conditions: 


Subroutine FORDIC loops over the entire grid to 
initialize each word of the cell dictionary. 

None 


Messages to Printer: 


None 






Element Type: FORTRAN Subroutine 


Element Name: FUNCT 

Entry Points: FUNCT 

Called By: GRAPHS 

External References: FXVSFY, VECTOR, ISOMET, MASSES, CONTUR, 

PLTSTR 

Common Blocks: SCCNTL 

Calling Sequence: (L) 

Processor: EDIT 


Variable 

Description 

L 

Pointer to type of plot. 


L = 2 

Velocity vectors 


L = 5 

Contours 


L = 6 

Principal stresses 


c*. 

II 

►4 

Maximum shear 


Other values for L are not valid. 


Element Description: Subroutine FUNCT selects the type of plot to be done 

on the SC 4020. 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: 


FXVSFY 


Element Type: FORTRAN Subroutine 


Entry Points: FXVSFY 

Called By: FUNCT 

External References: ABORT 

Common Blocks: None 

Calling Sequence: None 

Processor: EDIT 


Variable 


Description 


Element Description: FXVSFY is a dummy routine. It was intended that 

FXVSFY eventually would become a routine which 
would perform parameter versus parameter plots. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: 

GEDATE 

Element Type: FORTRAN Subroutine 

Entry Points: 

GEDATE 


Called By: 

FILE 



External References: None 

Common Blocks: None 

Calling Sequence: DATE 


Processor: 


I.C. , INTEG, EDIT 


Variable 

Description 

DATE 

Two words which when written in 2A6 format 
supply the current date. 


Element Description: 

Abort Conditions: 
Messages to Printer: 


Now a dummy routine, GEDATE was originally a 
SLEUTH routine written for EXEC II . 

None 

None 
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Element Name: GENBND Element Type: FORTRAN Subroutine 

Entry Points: GENBND 

Called By: MAIN1 

External References: ABORT, MASSIJ, PACUNP, PRTBPT 

Common Blocks: SUBSBK, BOUNDY, CONSTS, CUTBLK, RC DATA , 

SHCOM, IOCNTL 

Calling Sequence: None 

Processor: None 



Element Description: GENBND reads the namelist BDSTRG to define new 

string points. Subroutine MASSIJ is called to 
determine the subscripts of the lattice and regular 
cells each string point is in. These subscripts 
are packed into the control word for the string 
point via subroutine PACUNP. Finally PRTBPT is 
called to print the string points ; - 

Abort Conditions: GENNIM - The user indicated that strain hardening 

was to be performed but no strain hardening string 
points were specified. 

GENBND - The strain hardening string points did 
not come first or there are no string points . 

Messages to Printer: GENBND prints out the namelist BDSTRG. There is 

a message associated with the first GENBND abort. 
If a string point is found to lie outside the limits 
of the problem a message indicating this is printed 
along with the values of the relevant data . 
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Element Name: GETCON 


Element Type: FORTRAN Integer 

Function Subroutine 


Entry Points: GETCON 

Called By: ARCEL, BDINIT , BDPOST, EXIST, FILLER, MAPMIX, 

PRNTCT, SETMAT, SETSEG, SET SID, STARVL . VOLCEL, 
CONTUR, LNEDIT, CBNDRY , SEGMNT , MASFLO 


External References: 
Common Blocks: 
Calling Sequence: 
Processor: 


SHFTR, TRACE, ABORT 
CONSTS, CBLOCK 
L,M 

I.C., EDIT 


Variable 

Description 

L 

M 

Subscript of cell dictionary word desired. 
Pointer to the type of cell information required. 


Element Description: 

Abort Conditions: 
Messages to Printer: 


Subroutine GETCON extracts items from the cell 
dictionary word. 

None 

N one 
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Element Name: GETIJ 

Element Type: FORTRAN Subroutine 

Entry Points: GETIJ 


Called By: ARECEL, VOLCEL, VARIBL 

External References: 

None 

Common Blocks: 

SUBSBK 

Calling Sequence: 

L, K 

Processor: 

I.C., INTEG, EDIT 



Variable 


Element Description: Given the subscript for the dictionary word of a cell, 

GETIJ returns the row and column indices for the cell 
through the common block SUBSBK. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: GETIME Element Type: FORTRAN Subroutine 

Entry Points: GETIME * 

Called By: FILE 


External References: 

None 

Common Blocks: 

None 

Calling Sequence: 

HOUR 

Processor: 

I.C. , INTEG, EDIT 


Variable 

Description 

HOUR 

Two words which when written in 2A6 format 
supply the time of day. 


Element Description: Now a dummy routine, GETIME was originally a 

SLEUTH routine written for EXEC II. 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: GRAPHS 


Element Type: FORTRAN Subroutine 


Entry Points: GRAPHS 

Called By: PLOTTR 

External References: PAPER, FUNCT 

Common Blocks: IOCNTL, SCCNTL 

Calling Sequence: None 

Processor: EDIT 



Element Description: Subroutine GRAPHS processes all SC4020 graphs 

for a drum file. GRAPHS calls PAPER to create the 
background of the plot followed by FUNCT to 
provide the plot foreground. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: GRIDVT 


Element Type: FORTRAN Subroutine 


Entry Points: GRIDVT 

Called By: MAIN1, MAIN 2, PLOTTR 

External References: SETSUB, IBET, ABORT 


Common Blocks: 
Calling Sequence: 
Processor: 


CON STS, SUBSBK, RCDATA, NETBLK 
None 

I.C., INTEG, EDIT 


Variable 


Description 


Element Description: Subroutine GRIDVT computes all the elements of the 

grid vector from the information in NETBLK. 


Abort Conditions: None of the abort conditions currently provided for in 

GRIDVT can occur. 

Messages to Printer: None 
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Element Name: HYPOEL Element Type: FORTRAN Subroutine 

Entry Points: HYPOEL 

Called By: COMP1 

External References: EXIST, SETSUB, DIVU , VOLCEL 

Common Blocks: CONSTS, MATBLK, DELSBK, IOCNTL, SEGBLK, 

SUBSBK, GRDATA, RC DATA , CUTBLK, DIVUBK 

Calling Sequence: None 


Processor: INTEG 


Variable 


Description 


l 


Element Description: Subroutine HYPOEL calculates the initial estimates 

of the deviatoric components of stress. The 
formulas used appear as equations 3.26 and 3.27 
in Reference 1 . The initial components of deviatoric 
stress are computed from the total components and 
an estimate of the pressure obtained from the trace 
of the o^j matrix at time N. If the stress cell does 
not exist both the estimated deviatoric components 
and the total stresses at time N are set to zero. 
This routine is not effective if both TWOMU and 
YIELD are zero. 

Abort Conditions: None 


Messages to Printer: None 
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Element Name: IBET 


Entry Points: IBET 

Called By: GRIDVT 


External References: 
Common Blocks: 
Calling Sequence: 
Processor: 


None 

None 

A, X, B, M 
I.C., EDIT 


Element Type: FORTRAN Integer 

Function Subroutine 


Variable 

Description 

A 

Left point of interval. 

B 

Right point of interval. 

X 

Point under consideration. 

M 

Dummy variable . 


Element Description: 


IBET determines if the integer X is between the 
the integers A and B. 


IBET = - 1 
= 0 
= 1 
= 2 
= 3 


X < A or X > B 
X >A and X < B 
X = A 
X = B 
X = A = B 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: INTPSH 


Element Type: FORTRAN Subroutine 


Entry Points: INTPSH 

Called By: MAIN1 , BDINIT, MAIN2 , MAIN 3 


External References: 
Common Blocks: 
Calling Sequence: 
Processor: 


PUTGET 

BOUNDY 

None 

I.C INTEG, EDIT 


Varia ble 


Description 


Element Description: Subroutine INTPSH initializes the push down list. 

This releases the entire push down list for use by 
the boundary processor. 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: INTRST Element Type: FORTRAN Subroutine 

Entry Points: INTRST 

Called By: BDNORM, BDSTRS 

External References: LINES, ABORT 

Common Blocks: None 

Calling Sequence: X, Y, IX, JX, XT, YT, IN, JN, ISIDE, XS, YS 

Processor: IC, INTEG 


Variable 

Description 

X, Y 

Coordinates of the string section to intersect. 

HH 

X 

£ 

Row and column indices for cell. 

XT, YT 

Coordinates of intersection . 

IN, JN 

I and J offsets for next cell. 

ISIDE 

Pointer to side of cell intersected. 

XS, YS 

Arrays containing cell sides. 


Element Description: INTRST calculates the intersection of a segment of 

the string and the sides of a cell. Should there be 
more than one intersection, the last one found is 
returned. 

Abort Conditions: INTRST - when no intersection is found. 


Messages to Printer: None 
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Element Name: 


ISOMET 


Element Type: FORTRAN Subroutine 


Entry Points: ISOMET 

Called By: FUNCT 

External References: ABORT 

Common Blocks: None 

Calling Sequence: None 

Processor: EDIT 


1 Variable 


Description 


Element Description: ISOMET is a dummy routine. It provides a place to 

add isometric plotting capabilities should this be 
required . 


Abort Conditions: ISOMET - if this subroutine is called. 

Messages to Printer: None 
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Element Name: LEDPLT _ Element Type: FORTRAN Subroutine, 

Entry Points: LEDPLT 

Called By: RESULT 

External References: PAPER 

Common Blocks: SCCNTL, IOCNTL, SCDEPS 

Calling Sequence: None 

Processor: EDIT 

Variable Description 


Element Description: LEDPLT plots selected cell variables for the linear 

editor. 

Abort Conditions: None 

Messages to Printer: When LEDPLT produces a plot it prints the message 

"GRAPH" followed by the plot number. 
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Element Name: LINES Element Type: FORTRAN Subroutine 

Entry Points: LINES 

Called By: FOLTES, INTRST 

External References: BET 


Common Blocks: 
Calling Sequence: 
Processor: 


None 

$, XI, Yl, X2 , Y2 , X, Y 
I.C., INTEG 


Variable 

Description 

$ 

Nonstandard return. 

XI 

X coordinate for first segment. 

Yl 

Y coordinate for first segment. 

X2 

X coordinate for second segment. 

Y2 

Y coordinate for second segment. 

X 

X coordinate of intersection 

Y 

Y coordinate of intersection 


Element Description: LINES calculates the intersection between two line 

segments. If the two segments intersecting the 
coordinates of the intersection are returned. If 
they do not intersect the nonstandard return is 
taken. 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: 
Entry Points: 
Called By: 

External Reference 
Common Blocks: 
Calling Sequence: 
Processor: 


LINPLT Element Type: FORTRAN Subroutine 

LINPLT 

CONTUR, LEDPLT, PAPER, PAPER 1 , PLTBND , PLTSTR, VECTOR 
s: None 

None 

XII, YI1 , XI2 , YI2 , ID1 , INCTR 
EDIT 


Variable 

Description 

XII, YI1 

x, y coordinate of start of line. 

XI2, YI2 

x, y coordinate of end of line. 

ID1 

Pointer to edge of plot exceeded by first point. 

INCTR 

Pointer to edge of plot exceeded by second point. 


Element Description: LINPLT draws lines for SC4020 plots. It first checks 

if either point is off the field of view. If so, the line 
is drawn up to edge of plot. 


Abort conditions: None 


Messages to Printer: None 
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Element Name: LNEDIT Element Type: FORTRAN Subroutine 

Entry Points: LNEDIT 

Called By: PLOTTR 

External References: GETCON, SOPEN, SREAD, SETSUB, BET, RESULT 

Common Blocks: SCCNTL, SUBSBK, CONSTS, RC DATA , SCDEPS, 

IOCNTL 

Calling Sequence: PLOT 

Processor: EDIT 



Element Description: LNEDIT was written for a specialized series of 

calculations involving the impact of a sphere on a 
thin plate. It extrapolates the current mass positions 
and velocities to calculate information about the 
scatter of material at a later time. 


Abort Conditions: None 


Messages to Printer: When a cell does not lie in the doner region selected 

or when the extrapolated values do not lie in the 
catacher region, diagnostic information is printed out. 
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Element Name: MAIN1 


Element Type: FORTRAN Main 
Routine 


Entry Points: None 

Called By: None 

External References: 

Common Blocks: 

Calling Sequence: 
Processor: 


LOADRG , COLD, DUMSUB, INTPSH, GRIDVT , GENBND, 
MOVER, DMSUB1 

SIOBLK 

None 

I.C. 


Varia ble 


Description 


Element Description: 
Abort Conditions: 


MAIN 1 is the main routine for the initial 
conditions . 

None 


Messages to Printer: 


None 






J 


Element Name: MAIN2 Element Type: FORTRAN Main 

-------- Routine 

Entry Points: None 

Called By: J None 

External References: COLD, CINOUT, GRIDVT, DATA 2 , INTPSH, 

COMPUT, PRNTCT, PRTBPT 

Common Blocks: SIOBLK, CONSTS, FRSTBK 

Calling Sequence: None 

Processor: INTEG 



Element Description: MAIN 2 is the main program for the integration processor. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: MAIN3 Element Type: FORTRAN Main 

Routine 

Entry Points: None 

Called By: None 

External References: COLD, DATA, INTPSH, PLOTTR 

Common Blocks: IOCNTL, CON STS, SENSWT 

Calling Sequence: None 

Processor: EDIT 

Variable Description 


Element Description: MAIN 3 is the main routine for the edit processor. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: MAPMIX 

Entry Points: MAPMIX 

Called By: PRNTIJ 

External References: GETCON, TITLE, SETSUB 

Common Blocks: RCTEMP, CONSTS, SUBSBK, IOCNTL 

Calling Sequence: None 

Processor: I.C., INTEG, EDIT 


Element Type: FORTRAN Subroutine 


Variable 


Description 



Element Description: Subroutine MAPMIX creates a display on the printer 

indicating if a cell is empty, cut, or full. Each 
print position represents a cell, a blank indicating 
a void, a 1 indicating a partially full cell and a 2 
indicating a completely full cell. 


Abort Conditions: None 

Messages to Printer: As Indicated 
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Element Name: MASFLO 

Entry Points: MASFLO 

Called By: COMP2 


Element Type: FORTRAN Subroutine 


External References: 
Common Blocks: 


Calling Sequence: 
Processor: 


GETCON, EXIST, SETSUB, SEGMNT, VOLCEL 

CONSTS, DELSBK, IOCNTL, SEGBLK, SUBSBK, 
C BLOCK, GRDATA , RCDATA, CUTBLK, RCTEMP 

IMASFL 

INTEG 


Variable 

Description 

IMASFL 

Column currently being computed. 


Element Description: Subroutine MASFLO calculates the flow terms for 

mass, momentum and energy. MASFLO then 
repartitions the completed cells. The mathematics 
used in MASFLO can be derived from the equations 
appearing in Section 3 . 2 . 2 . 1 of Reference 1 . 


Abort Conditions: None 

Messages to Printer: MASFLO prints out cycle by cycle totals of mass, 

momentum, and energy. If during a calculation 
negative energy was set to zero, the total amount 
of negative energy for the cycle is printed following 
the message "TOTAL NEGATIVE ENERGY SUPPRESSED" . 
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Element Name: 
Entry Points: 
Called By: 

External Reference 
Common Blocks: 
Calling Sequence: 
Processor: 


MASSES 
MASSES 
FUNCT 
s: ABORT 

None 
None 
EDIT 


Element Type: FORTRAN Subroutine 


Variable 


Description 


Element Description: MASSES is currently a dummy subroutine. It was 

intended that additional plotting capabilities would 
be inserted here eventually. 


Abort Conditions: MASSES - if subroutine is called. 

Messages to Printer: None 


72 






Element Name: MASSIJ Element Type: FORTRAN Subroutine 

Entry Points: MASSIJ 

Called By: GENBND, BDMOVE, COLPTS 

External References: SETSUB 

Common Blocks: SUBSBK, CON STS, RC DATA 


Calling Sequence: 
Processor: 


$, K, XM, YM 
I.C., INTEG 


Variable 

Description 

$ 

Nonstandard return. 

K 

= 0 All four subscripts are obtained. 


= 1 Only I and J are obtained. 


= 2 Only II and LJ are obtained. 

XM 

X coordinate of string point . 

YM 

Y coordinate of string point. 


Element Description: Subroutine MASSIJ determines the subscripts of the 

regular and lattice cell from the mass portion 
for a particular string point. 


Abort Conditions: None 

Messages to Printer: None 
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Element Type: FORTRAN Subroutine 


Element Name: MOVER 

Entry Points: MOVER 

Called By: MAIN 1 , COMPUT 

External References: BDINIT, BONDPR, FORDIC , BDPOST 

Common Blocks: BOUNDY 

Calling Sequence: None 

Processor: I.C., INTEG 



Element Description: Subroutine MOVER processes the boundary string 

after the string points have been moved. It is the 
driver routine for the string processor used merely 
to call BDINIT, BONDPR, FORDIC and BDPOST. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: 
Entry Points: 
Called By: 

External Reference 
Common Blocks: 
Calling Sequence: 


OUTPUT Element Type: 

OUTPUT 

DMSUB1 

s: PRNTCT, PRTBPT, RESTOU 

None 
None 


Processor: 


I.C. 


FORTRAN Subroutine 


Variable 


Description 


Element Description: Subroutine OUTPUT is used by initial conditions to 

produce printed output of the control dictionary, 
string positions, and the cell variables. It also 
writes the initial conditions restart. 

Abort Conditions: None 

Messages to Printer: None 


75 





Element Name: PACKER Element Type: FORTRAN Subroutine 

Entry Points: PACUNP 

Called By: BONDPR, GENBND, PRTBPT, COLPTS 

External References: SHIFTL, SHIFTR 


Common Blocks: 
Calling Sequence: 
Processor: 


None 

T, S, MODE 
I.C., INTEG 


Variable 

Description 

MODE 

= 1; Packing items in array "s” into word "T". 

= 2; Unpacking items from word "T" into array "s 11 . 


Element Description: Subroutine PACKER packs and unpacks the subscripts 

from the dictionary word for the string points . 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: PAPER 


Element Type: FORTRAN Subroutine 


Entry Points: PAPER 

Called By: GRAPHS , LEDPLT 

External References: None 

Common Blocks: SCDEPS, CONSTS, SCCNTL 

Calling Sequence: I 

Processor: EDIT 



Element Description: Subroutine PAPER draws the background for an SC 4 020 

plot. PAPER first advances the frame. It then provides 
the standard titles, labels, and grid. If the plot is 
a contour plot, additional labeling is provided to 
distinguish the different contours. If the plot is a 
velocity or stress field, a characteristic vector length 
is plotted. 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: PHASEA Element Type: FORTRAN Subroutine 

Entry Points: PHASEA 

Called By: COMP1 

External References: CBN DRY, SETSUB, VOLCEL, SETFOR, COLPTS 

Common Blocks: CONSTS, DELSBK, IOCNTL, SEGBLK, SUBSBK, 

C BLOCK,. GRDATA , RC DATA , CUTBLK, RCTEMP, 

TOLBLK 

Calling Sequence: IPHSA 

Processor: INTEG 


Variable 

Description 

IPHSA 

Index to column being computed. 


Element Description: Subroutine PHASEA is called during the initial part 

of the integration computations to calculate the 
energy and velocity estimates. PHASEA calls 
SETFOR to determine the forces on the sides of the 
cell. The equations for calculating the estimated 
velocities and energies appear as equations 3.20, 
3.21, and 3.25 in Reference 1. Incidental to these 
calculations PHASEA calls COLPTS to move the 
string points in a column. 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: PHI Element Type: FORTRAN Real 

Function Subroutine 

Entry Points: PHI 

Called By: PLASFL 

External References: None 

Common Blocks: None 

Calling Sequence: None 

Processor: INTEG 



Element Description: Subroutine PHI is a special function used in 

conjunction with the Perzyna model. A description 
of the Perzyna model can be found in Section 2.2 of 
Reference 2 . 


Abort Conditions: None 


Messages to Printer: None 
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Element Name: 
Entry Points: 
Called By: 

External Reference 
Common Blocks: 

Calling Sequence: 
Processor: 


PLASFL Element Type: FORTRAN Subroutine 

PLASFL 

vise OS 

s: ABORT, YLD, FCUT, THETA, PHI 

MATBLK, CONSTS, RCDATA, TOLBLK, BOUNDY, SHCOM, 
CUTBLK, SUBSBK, DIVUBK 

SXX, SYY, STT, TXY 

INTEG 


Varia ble 

Description 

SXX 

SYY 

STT 

The deviatoric stresses. 

TXY 

Element Description: 

Subroutine PLASFL imposes the yield conditions 
upon the deviatoric stresses and calculates the 
incremental change in the strain rate. Currently 
either the von Mises or Perzyna models may be 
selected. 

Abort Conditions: 

PLASFL if there are more than 50 string points 
in a particular stress cell. 

Messages to Printer: 

If the Perzyna model is selected there are three 
messages which may be printed in conjunction 
with the Newtons' method. These messages 
start off 


"XN TOO LARGE" 


"NEWTONS METHOD FAILS TO CONVERGE" 
"F FROM NEWTON IS TOO LARGE" 
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Element Name: PLOTTR Element Type: FORTRAN Subroutine 

Entry Points: PLOTTR 

Called By: MAIN 3 

External References: ABORT, FILE, GRIDVT, LNEDIT, GRAPHS -S 


Common Blocks: 
Calling Sequence: 
Processor: 


IOCNTL, SCCNTL, SCDEPS 

N one 

EDIT 


Variable 


Description 



Element Description: Subroutine PLOTTR is one of the main routines in the 

in the SC4020 processor. PLOTTR reads the data 
which controls both the plots and the restarting . 
Subroutine GRAPHS is called to actually produce 
the plots . 


Abort Conditions: PLTR.l - if the number of files to edit is zero and 

DOALL = .FALSE. 

PLTR.2 - NGRAPH .LE.O 


Messages to Printer: PLOTTR lists the input data. 
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Element Name: PLTBND 


Element Type: FORTRAN Subroutine 


Entry Points: PLTBND 

Called By: CONTUR, PLTSTR, VECTOR 

External References: 

Common Blocks: BOUNDY, SCCNTL 

Calling Sequence: M 

Processor: EDIT 


Variable 

Description 

M 

Type of graph being plotted. 


Element Description: Subroutine PLTBND plots string points. The points 

which are associated with real strings are connected 
with straight lines. 


Abort Conditions: None 


Messages to Printer: None 


82 




Element Type: FORTRAN Subroutine 


Element Name: PLTSTR 

Entry Points: PLTSTR 

Called By: FUNCT 

External References: EXIST, SOPEN, SREAD, SETSUB 

Common Blocks: CONSTS, SUBSBK, RC DATA , GRDATA , SCCNTL 

Calling Sequence: L, MM 

Processor: EDIT 


Variable 

Description 

L 

Graph Number 

MM 

= 1; Principal stresses are plotted. 

= 2; Maximum shear stresses are plotted. 

i 

1 - 


Element Description: PLTSTR plots either the principal stresses or the 

maximum shear. PLTSTR makes one pass over the 
drum file . 


Abort Conditions: None 


Messages to Printer: None 
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Element Name: 
Entry Points: 
Called By: 

External Reference 
Common Blocks: 
Calling Sequence: 
Processor: 


PNTSET 
PNTSET 
COMP1 
s: PACUNP 

BOUNDY 
None 
INTEG 


Element Type: FORTRAN Subroutine 


Variable Description 


Element Description: PNTSET sets each string control word to indicate 

that the string point has not been moved. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: 
Entry Points: 
Called By: 

External Reference 
Common Blocks: 
Calling Sequence: 
Processor: 


PRESS Element Type: 

PRESS 

STRESS 

s: TBFT02 

PCNBLK, EOSBLK, MATBLK 
IS 

I.C., INTEG 


FORTRAN Subroutine 


Variable 

Description 

IS 

Material number being considered. 


Element Description: Subroutine PRESS is called by subroutine STRESS to 

determine the pressure at a lattice point. Currently 
it assumes that the equation of state TBFT02 is 
being used. The density and specific internal energy 
are inputted through the common block EOSBLK as 
variables labeled ETA and EN respectively. The 
pressure and sound speed squared are in turn out- 
putted in EOSBLK as the variables PRES and VL. 


Abort Conditions: None 


Messages to Printer: None 
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Element Name: PRESTI Element Type: FORTRAN Subroutine 

Entry Points: PRESTI 

Called By: COMPUT 

External References: None 

Common Blocks: None 

Calling Sequence: None 

Processor: INTEG 


Variable Description 


Element Description: PRESTI supplies externally applied pressure to subroutine 

COMPUT. PRESTI currently has no pressure function. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: PRNTCT 


Element Type: FORTRAN Subroutine 


Entry Points: PRNTCT 

Called By: OUTPUT, MAIN 2 

External References: GETCON, SETSUB, TITLE 

Common Blocks: CONSTS, IOCNTL, SUBSBK, CBLOCK 

Calling Sequence: None 

Processor: I.C., INTEG 



Element Description: Subroutine PRNTCT prints the contents of the cell 

dictionary. Lines that are equal by columns are 
not printed . 


Abort Conditions: None 


Messages to Printer: None 
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Element Name: PRNTIJ Element Type: FORTRAN Subroutine 

Entry Points: PRNTIJ 

Called By: FILE 

External References: EXIST, SETSUB, TITLE, ARECEL, VOLCEL, MAPMIX 

Common Blocks: CONSTS, DELSBK, IOCNTL, SUBSBK, SEGBLK, 

C BLOCK, GRDATA , CUTBLK, RCDATA, 

RCTEMP 

Calling Sequence: IPRT 

Processor: I.C., INTEG, EDIT 


Variable 

Description 

IPRT 

Current column number being printed. 


Element Description: Subroutine PRNTIJ prints the values of regular cell 

variables for a given column. This is followed by 
the values of the stress cell variables for that 
column. PRNTIJ is used as a subroutine to RESTOU 
to avoid making multiple passes over the drum. 
PRNTIJ also totals the momentums, kinetic energy, 
internal energy, and the total energy. These totals 
are printed out following the stress cell variables 
for the last column. 

For the integration or COMPUT processor PRNTIJ, 
checks on activities before printing. A regular cell 
is considered active if either component of velocities 
or the internal energy for that cell is non-zero. A 
stress cell is considered active if any one of the 
oij are non-zero. 

Abort Conditions: None 


Messages to Printer: As indicated. 
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Element Name: 
Entry Points: 
Called By: 

External Reference 
Common Blocks: 
Calling Sequence: 
Processor: 


PRTBPT Element Type: FORTRAN Subroutine 

PRTBPT 

GENBND, OUTPUT, MAIN2 
s: TITLE, PACUNP 

IOCNTL, BOUNDY, SHCOM, CUTBLK 
None 

I.C., INTEG 


Variable 


Description 


Element Description: PRTBPT prints the positions of all the string points. 

Abort Conditions: None 


Messages to Printer: 


As indicated. 






Element Name: PUTCON Element Type: FORTRAN Subroutine 

Entry Points: PUTCON 

Called By: BDINIT, FILLER, FORDIC , SETMAT, SETSEG , SETSID, 

STARVL 

SHFTL 

CONSTS, C BLOCK 
A, L, M 
I.C . , INTEG 


Varia ble 

Description 

A 

The data item to enter into the dictionary. 

L 

Location of the dictionary word. 

M 

= 0; Entire word "A" is set into the dictionary. 

- 

/ 0; Then M is the key for the item to be placed 
in the dictionary. 


Element Description: PUTCON puts control information for a regular cell or 

lattice cell into the cell dictionary. PUTCON is the 
counter part to GETCON. 

Abort Conditions: None 

Messages to Printer: None 


External References: 
Common Blocks: 
Calling Sequence: 
Processor: 
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Element Name: PUTGET 

Entry Points: PUTGET 

Called By:. BDINIT, INTPSH/ SETSEG, STARVL 

External References: ABORT 

Common Blocks: CUTBLK, BOUNDY 

Galling Sequence: LOC , MODE 

Processor: I.C., INTEG 


Element Type: FORTRAN Subroutine 


Variable 

Description 

MODE 

= 1; "LOC" is obtained and set by this routine. 

= 2; "LOC" is released and next available 

location is set to the value of "LOC". 


Element Description: PUTGET manipulates the pushdown list in CUTBLK. 


Abort Conditions: GET - There is an attempt to get a location in the 

pushdown lists and then are no more available. 

PUT - There is an attempt to release a location whose 
value is 0 or greater than the number of available 
locations. 


Messages to Printer: None 
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Element Name: REGION Element Type: FORTRAN Subroutine 

Entry Points: REGION 

Called By: DMSUB1 

External References: SOPEN , SREAD, SETSUB, VOLCEL, SRITER, 

SCLOSE 

Common Blocks: CONSTS, SUBSBK, IOCNTL, GRDATA, RCDATA, 

REGBLK, MATBLK, EOSBLK 

Calling Sequence: MODE 


Processor: 


I.C. 


Variable 

Description 

MODE 

Dummy and argument. 


Element Description: REGION sets the initial values of the cell variables. 

REGION leaves the machine refreshed with a new file 
on drum. -• 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: RESTIN Element Type: FORTRAN Subroutine 

Entry Points: RESTIN 

Called By: CINOUT 

External References: FILE 

Common Blocks: IOCNTL 

Calling Sequence: None 

Processor: INTEG 



Element Description: RESTIN controls the refreshing process from the restart 

tape. As its input RESTIN reads the namelist RSTART. 
Subroutine FILE is called both to position the tape and 
to read the desired file. 

Abort Conditions: None 

Messages to Printer: The namelist RSTART is printed. In addition 

information identifying the files being read by 
subroutine FILE is provided. 
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Element Name: RESTOU 

Entry Points: RESTOU 

Called By: OUTPUT, CINOUT 

External References: FILE 

Common Blocks: IOCNTL 


Calling Sequence: None 

Processor: I.C., INTEG 


Element Type: FORTRAN Subroutine 


Variable Description 


Element Description: RESTOU writes a restart file on the output tape. 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: RESULT 


Element Type: FORTRAN Subroutine 


Entry Points: RESULT 

Called By: LNEDIT 

External References: LEDPLT 

Common Blocks: CON STS, SCCNTL, GRDATA , RCDATA, SCDEPS, 

IOCNTL 

Calling Sequence: IGOTO, PLOT 

Processor: EDIT 


Varia ble 

Description 


= 1; Additional printout is taken 

IGOTO 

^ = 2; Only key printout is taken 

PLOT 

= .TRUE.; Subroutine LEDPLT, forcing plotting 
of linear editor. 


Element Description: Result is the control routine for print and plot in the 

linear editor. 


Abort Conditions: None 


Messages to Printer: RESULT prints out the results of the linear 

extrapolation. This is the mass, momenta, and 
energies in each of the catcher regions. 
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Element Name: SAVSUB Element Type: FORTRAN Subroutine 

Entry Points: SAVSUB 

Called By: ARECEL, VOLCEL, VARIBL 

External References: None 

Common Blocks: SUBSBK 

Calling Sequence: MODE 

Processor: I.C., INTEG, EDIT 


Variable 

Description 

MODE 

= 1; Items in SUBSBK are stored in local 
storage OU. 

= 2; Items in SUBSBK are restored in local 
storage IN . 


Element Description: SAVSUB saves and restores subscript information 

in common block SUBSBK. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: SEGMNT 


Element Type: FORTRAN Subroutine 


Entry Points: SEGMNT 

Called By: MASFLO, SETDEL, SETFOR 

External References: GETCON 

Common Blocks: CONSTS, IOGNTL, SEGBLK, SUBSBK, RC DATA , 

CUTBLK 

Calling Sequence: (L, M) 

Processor: INTEG 


Variable 

Description 

L 

= 1; Horizontal segment. 


= 2; Vertical segment. 

M 

= 1; Regular cell. 


= 2; Lattice cell. 


Element Description: Subroutine SEGMNT determines the type of segment 

and the end points of the segments for either the 
vertical or horizontal side of a stress or regular 
cell. The output from this routine appears in common 
block SEGBLK. The output from SEGMNT is as follows: 

NOSEG - The numbers of segments containing mass. 
NOCUT - The number of times a cell is cut. 



XA 

XB 

- X coordinates 

- Y coordinates 

of the cuts . 
of the cuts . 

Abort Conditions: 

None 



Messages to Printer: 

None 
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Element Name: 


SELDT 


Element Type: FORTRAN Subroutine 


Entry Points: 
Called By: 

External Reference 
Common Blocks: 
Calling Sequence: 
Processor: 


SELDT 

COM PUT 

s: BET, ABORT 

CONSTS, IOCNTL 

None 

INTEG 


Variable 


Description 


Element Description: SELDT monitors the stability conditions for a STEEP 

solution and adjusts the integration time step DT 
to hold the solution stable. The inputs to SELDT 
are CMAX and VMAX. CMAC is calculated by sub- 
routine STRESS. VMAX is calculated by subroutine 
MASFLO. The quantity SMAX is set to the maximum 
of CMAX and VMAX. DT is then decreased if SMAX 
is greater than SI or increased if it is less than S2. 

The following formulas define CMAX and VMAX: 

2 

CMAX = [ ( part:i - c l e Velocity + Dilatational Wave Speed) * At -j 


2 

VMAX = [ Maximum Particle Velocity * At ^ 
^Minimum (AX, AY) 


Abort Conditions: None 


Messages to Printer: None 
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Element Name: 
Entry Points: 
Called By: 


SETDEL 

SETDEL 

STRFLO 


Element Type: FORTRAN Subroutine 


External References: 
Common Blocks: 

Calling Sequence: 
Processor: 


EXIST, SETSUB, SEGMNT 

CONSTS, DELSBK, IOCNTL, SUBSBK, SEGBLK, 
C BLOCK, GRDATA , RC DATA , RCTEMP , CUTBLK 

II, JJ, MODE 

INTEG 


Variable 

Description 

II 

Column pointer for lattice cell. 

JJ 

Row pointer for lattice cell. 


= 1; Stress flow across horizontal surface being 
ca lculated . 

MODE 

= 2; Stress flow across vertical surface being 
calculated. 


= 3; Stress flow across both surfaces being 
calcula ted. 


Element Description: SETDEL computes the stress flow across the 

selected surfaces of a lattice cell. The actual 
flow for each of the components is returned through 
common block DELSBK. 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: 
Entry Points: 
Called By: 

External Reference 

Common Blocks: 

Calling Sequence: 
Processor: 


SETFOR Element Type: FORTRAN Subroutine 

SETFOR 
P HA SEA 

s: EXIST, SURPRE, SETSUB, SEGMNT, BOUNST, 

ARECEL 

CONSTS, DELSBK, IOCNTL, SUBSBK, SEGBLK, 

GRDATA , CUTBLK, RCDATA, RCTEM? 

II, JJ, MODE 

INTEG 


Variable 

Description 

II 

Column pointer for regular cell. 

JJ 

Row pointer for regular cell. 


= 1; Forces for horizontal surface are set. 

MODE 

= 2; Forces for vertical surface are set. 


= 3; Forces for both surfaces and force into 
plane are set. 


Element Description: SETFOR sets the horizontal and vertical forces on the 

surfaces of a regular cell. The final forces are 
returned through RCTEMP. 

Abort Conditions: None 


Messages to Printer: None 
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Element Name: SETMAT 

Entry Points: SETMAT 

Called By: BDPOST 

External References: GETCON, SETSUB, ABORT, PUTCON 

Common Blocks: CONSTS, SUBSBK, CUTBLK 

Calling Sequence: None 

Processor: I.C., INTEG 


Element Type: FORTRAN Subroutine 


Variable 


Description 


Element Description: SETMAT is used during the boundary processing to 

set the cell dictionary word. 

Abort Conditions: SETMAT - If a negative value was extracted from 

the cell dictionary. 

Messages to Printer: None 
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Element Name: SETSEG 


Element Type: FORTRAN Subroutine 


Entry Points: SETSEG 

Called By: BDNORM, BDSTRS 

External References: GETCON, SETSUB, PUTGET, PUTCON, ABORT 


Common Blocks: 
Calling Sequence: 
Processor: 


CON STS, SUBSBK, CUTBLK 
IFLG, I, J, X, Y, SIDE, M 
I.C. , INTEG 


Variable 

Description 

I, J 

Column and row indices of cell under consideration. 

X, Y 

Intersection points . 

SIDE 

Side which is cut. 

M 

= 1; Regular cell, = 2; Stress cell. 

IFLG 

Not used. 


Element Description: SETSEG allocates space for a cut cell and stores 

information in CUTBLK. 

Abort Conditions: SETSG2 - If a boundary point lies on a cell side. 

Messages to Printer: SETSEG prints three types of diagnostic messages. 

The first indicates that a string point was found to 
lie on a cell side . The other two indicate that a 
LOC was or was not released. 
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Element Type: 


Element Name: SETSID 

Entry Points: SETSID 

Called By: BDPOST 

External References: GETCON, ABORT, PUTCON 

Common Blocks: SCANBK 

Calling Sequence: None 

Processor: I.C., INTEG 


FORTRAN Subroutine 


Variable 


Description 


Element Description: SETSID scans over the segment types previously set 

in the cell dictionary and fills in particular key items. 
SETSID sweeps over the dictionary line by line in the 
direction specified by the calling routine through 
common block SCANBK. 


Abort Conditions: SETSID 

Messages to Printer: None 


Either an initial key was found to be 
negative or an illegal key was found. 
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Element Name: 


SETSUB 


Element Type: FORTRAN Subroutine 


Entry Points: SETSUB 

Called By: BDNORM, BDPOST, BDSTRS , COLD, DATA, FILE, FILLER, 

FORDIC, GRIDVT, MAPMIX, MASSIJ, PRNTIJ, PRNTCT, 
REGION, SETMAT, SETSEG, STARVL , BDMOVE , HYPOEL, 
PHASEA, SETDEL, SETFOR, STRESS, STRFLO , MASFLO, 
CONTUR, LNEDIT, PLTSTR, VARIBL, VECTOR 


External References: 
Common Blocks: 
Calling Sequence: 
Processor: 


ABORT 

CONSTS, SUBSBK 
K 

I.C., INTEG, EDIT 


Variable 

Description 

K 

= 0; Initialize offset values for problem. 

= 1; Cell reference IJ and dictionary reference IJD are 
set from cell indices I and J. 

= 2; Lattice cell information is from LI and LJ. 

= 3; Lattice cell information is from I and J. 

= 4; Cell and lattice cell information is from I and J. 

= 5; Cell and lattice cell information is from LI and LJ. 

= 6; I and J are set from LI and LJ. 

= 7; Regular cell information is set from LI and LJ. 


Ellement Description: Subroutine SETSUB serves two functions . First for 

K = 0 it initializes the quantities which allow the 
programmer to refer to offset cells in GRDATA and 
in the cell dictionary. For the other values of K 
items are set in common block SUBSBK. SETSUB 
and SIO are the only two routines which have direct 
information about the rotation of the buffers in 
GRDATA . 


Abort Conditions: SETSUB - If the column or row index causes a 

violation of the code dimensions. 


Messages to Printer: None 
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Element Name: 
Entry Points: 
Called By: 

External Reference 
Common Blocks: 
Calling Sequence: 
Processor: 


SETUP 
SETUP 
STARVL 
s: None 

CUTBLK 

L 

I.C., INTEG 


Element Type: FORTRAN Subroutine 


Variable 

Description 

L 

Superscript for variables in CUTBLK. 


Element Description: 

Abort Conditions: 
Messages to Printer: 


SETUP initializes values in common block CUTBLK 
subscripted by L. 

None 

None 
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Element Name: 

Entry Points: 
Called By: 

External Reference 
' Common Blocks: 
Calling Sequence: 
Processor: 


SHFTL 

SHFTL 
PUTCON 
s: None 

None 

IWORD, NBITS 
I.C. 


Element Type: FORTRAN Integer 

Function Subroutine 


Variable 

Description 

IWORD 

NBITS 

Word to be shifted. 

Number of bits to be shifted. 


Element Description: SHFTL currently outputs the contents of IWORD, 

shifted left NBITS, by multiplying by powers of 
two. Since this is inefficient, SHFTL should be 
replaced by a SLEUTH routine or the references 
to it should be changed to use an appropriate 
system routine if one is available. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: SHFTR 

Entry Points: SHFTR 

Called By: GETCON, PACKER 

External References: None 

Common Blocks: 


Element Type: FORTRAN Integer 

Function Subroutine 


Calling Sequence: 


None 

IWORD, NBITS 


Processor: 

I.C., 1NTEG, EDIT 

Variable 

Description 

IWORD 

Word being shifted. 

NBITS 

Number of bits to shift. 


Element Description: SHFTR shifts right. Analogous to SHFTL. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: SHIFTL Element Type: FORTRAN Integer 

Function Subroutine 

Entry Points: SHIFTL 

Called By: PACKER 

External References: None 

Common Blocks: None 

Calling Sequence: IWORD, NBITS 


Processor: 

I.C., INTEG 

Variable 

Description 

IWORD 

Word to be shifted. 

NBITS 

Number of bits to shift. 


Element Description: Same as SHFTL 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: SHIFTR 


Entry Points: SHIFTR 

Called By: PACKER 

External References: None 

Common Blocks: None 


Calling Sequence: IWORD, NBITS 

Processor: I.C., INTEG 


Element Type: FORTRAN Integer 

Function Subroutine 


Variable 

Description 

IWORD 

Word to be shifted. 

NBITS 

Number of bits to shift. 

Element Description: 

Same as SHFTR. 

Abort Conditions: 

None 

Messages to Printer: 

None 


109 



Element Name: SIO Element Type: FORTRAN Subroutine 

Entry Points: SOPEN, SREAD, SRITER, SCLOSE 

Called By: FILE, REGION, CONTUR, LNEDIT, PLTSTR, VECTOR, 

COMF1, COMP2 

External References: ABORT, NTRAN 

Common Blocks: GRDATA , IOCNTL, CON STS, SIOBLK 

Calling Sequence: None 

Processor: I.C., EDIT, INTEG 



Element Description: Subroutine SIO controls the reading and writing of 

the scratch drum files. To open a scratch file and 
initialize reading a call to SOPEN must be made. 

A call to SREAD insures that the previously initiated 
read is completed and issues a new read. A call 
'i to SWRITE waits for completion of the previous write 
and then initiates the next write. A call to CLOSE 
forces out all waiting writes, writes an end of file, 
and swaps the input and output units. 


Abort Conditions: NOFILE An attempt was made to open a file which 

was never created. 

EROPEN An attempt was made to read a file which 
was not opened. 

IOERR A read or write transmission error occurred. 


Messages to Printer: None 
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Element Name: STARVL Element Type: FORTRAN Subroutine 

Entry Points: STARVL 

Called By: BDNORM, BDSTRS 

External References: GETCON, SETSUB, PUTGET, PUTCON, SETUP 

Common Blocks: CONSTS, SUBSBK, CUTBLK 

Calling Sequence: I, J, A, V, MATL, M 


Processor: 

I.C. , INTEG 

Variable 

Description 

M 

= 1; For regular cells. = 2; For lattice cells. 

IJ 

Subscripts to the cell in question. 

A, V 

Incremental area and volume 

MATL 

= .TRUE. , Point LPT is set for lattice cells. 


Element Description: STARVL stores the area and volume for a cut cell. If 

the cell had a previous area and volume A and V 
... . are added in. Either CMAT or LMAT is set to 1 to 

indicate the cell is cut. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: STPOST 


Element Type: FORTRAN Subroutine 


Entry Points: STPOST 

Called By: COMPUT 

External References: None 

Common Blocks: SHCOM, CUTBLK 

Calling Sequence: None 

Processor: INTEG 



Element Description: STPOST updates the strain hardening at the string 

points which carry strain hardening. The change 
in strain hardening is set to zero in preparation for 
the next cycle . 

Abort Conditions: None 


Messages to Printer: None 
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Element Name: STRESS Element Type: FORTRAN Subroutine 

Entry Points: STRESS 

Called By: COMP2 

External References: EXIST, SETSUB, SURPRE, VOLCEL, ABORT, PRESS, 

VISCOS 

Common Blocks: DELSBK, DIVUBK, EOSBLK, MATBLK, CONSTS, 

IOCNTL, SUBSBK, SEGBLK, CBLOCK, GRDATA , 

D STARS, RC DATA , SHCOM, BOUNDY, TOLBLK 

i Calling Sequence: ISTRS 

Processor: INTEG 



Element Description: Subroutine STRESS is called to calculate the final 

stress. It is the last major routine in the integration 
process. The inputs to STRESS are the deviatoric 
components of stress for the column ISTRS. An 
effective energy and density are calculated from the 
four surrounding regular cells. The effective energy 
is calculated using mass weighting. The effective 
density is calculated as the sum of the masses of the 
surrounding regular cells divided by the sum of their 
volumes. The equation of state is then used to produce 
the pressure and sound speed of the material. VISCOS 
is called next to impose the yield condition and to 
calculate the artificial viscosities. The resulting 
pressure and viscosity are added to the constrained 
deviators to produce the final stresses. 

Abort Conditions: STRES2 A lattice point exists but none of the 

surrounding regular cells have mass. 

STRES3 A lattice point exists but none of the 

surrounding regular cells have volume. 

Messages to' Printer: PRESSURE TOO LARGE followed by identifying 

information indicates P ^ PMAX. 

MASS SUPPRESSED followed by identifying 
information indicates the effective densities is 
greater than an inputted maximum supplied by user. 
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Element Name: STRESS 


Messages to Printer: Continued 


PARTVL TOO LARGE followed by the value of the 
particle velocity indicates that the cell particle 
velocity was greater than the value of TOLNEW(3) 
supplied by the user. This cell will not be allowed 
to enter into the stability calculations. 


Following the calculations for the last column in the 
problem, the peak values, stresses, and pressure 
are printed for each cycle . 
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Element Name: 
Entry Points: 
Called By: 

External Reference 
Common Blocks: 

Calling Sequence: 
Processor: 


STRFLO Element Type: FORTRAN Subroutine 

STRFLO 

COMP2 

s: EXIST, SETSUB, SETDEL, VOLCEL 

CONSTS, DELSBK, IOCNTL, SUBSBK, SEGBLK, 

C BLOCK, GRDATA , RC DATA , CUTBLK, RCTEMP, 
MATBLK 

ISTRFL 

INTEG 


Variable 

Description 

ISTRFL 

Index of column to compute. 


Element Description: STRFLO is one of the major subroutines in the second 

part of the integration processor. STRFLO completes 
the flow terms in the deviatoric stress components 
for a given column and thus determines the final 
deviatoric components of stress. STRFLO is not 
executed for a purely hydrodynamic solution. 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: SURPRE 

Entry Points: SURPRE 

Called By: SETFOR, STRESS 

External References: None 

Common Blocks: None 

Calling Sequence: (I, J) 


Element Type: FORTRAN Real 

Function Subroutine 


Processor: 

INTEG 

Variable 

Description 

I 

Column pointer. 

J 

Row pointer. 


Element Description: SURPRE supplies surface pressure when required, 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: TBFT02 Element Type: FORTRAN Subroutine 

Entry Points: TBFT02 

Called By: PRESS 

External References: BET 

Common Blocks: EOSBLK, PCNBLK 

Calling Sequence: None 

Processor: I.C., INTEG 


Variable Description 



Element Description: TBFT02 is a tabular equation of state for 

aluminum . 

Abort Conditions: None 

Messages to Printer: None 


117 






Element Name: THETA - - Element Type: FORTRAN Real 

Function Subroutine 

Entry Points: THETA 

Called By: PLASFL 

External References: None 

Common Blocks: None 

Calling Sequence: EP 

Processor: INTEG 


Variable 

Description 

EP 

Plastic strain, e^. 


Element Description: THETA is a specialized function used in conjunction 

with the Perzyna work hardening model. 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: TITLE Element Type: FORTRAN Subroutine 

Entry Points: TITLE 

Called By: DATA, MAPMIX, PRNTIJ, PRNTCT, PRTBPT, DATAL 

External References: IOCNTL,. CONSTS 

Common Blocks: None 

Calling Sequence: None 

Processor: I.C., INTEG 



Element Description: TITLE writes a standardized title for the 

print file . 

Abort Conditions: None 

Messages to Printer: None 
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Element Name: VARIBL . Element Type: FORTRAN Real 

Function Subroutine . 

Entry Points: VARIBL 

Called By: CONTUR 

External References: EXIST, VOLC EL, SAVSUB, GETIJ, SETSUB 

Common Blocks: SUBSBK, CONSTS, GRDATA 

Calling Sequence: IJ, TYPE 

Processor: EDIT 


Variable 

Description 

IJ 

Is a subscript either to the GRDATA buffer or to the 


cell dictionary depending on the value for type. 


= 1; Specific internal energy. 


= 2; Density. 


= 3; Dynamic pressure. 

TYPE 

= 4; Kinetic energy. 


= 5; U momentum. 


= 6; V momentum. 


= 7; P at center of mass. 


Element Description: The function VARIBL returns the value of a regular cell 

variable to be used in the contour plots. 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: VECTOR Element Type: FORTRAN Subroutine 

Entry Points: VECTOR 

Called By: FUNCT 

External References: EXIST, SOPEN, PLTBND , SREAD , SETSUB, LINPLT 

Common Blocks: CONSTS, SUBSBK, GRDATA , SCCNTL 

Calling Sequence: L 

Processor: EDIT 


Variable 

Description 

L 

Number of the plot. 


Element Description: VECTOR plots the velocity vector field. VECTOR 

also calls PLTBND to plot the boundary. It 
should be noted that VECTOR requires equal X-Y 
scaling be specified in PAPER. VECTOR makes one 
pass over the drum file. 


Abort Conditions: None 

Messages to Printer: None 
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Element Name: 
Entry Points: 
Called By: 

External Reference 
Common Blocks: 

Calling Sequence: 
Processor:. 


VISC OS Element Type: FORTRAN Subroutine 

VISCOS 

STRESS 

s: DIVU , PLASFL 

EOSBLK, CON STS, SUBSBK, GRDATA, RC DATA , DIVUBK 
D STARS 

IT, JT, SXX, SYY, STT TXY 
INTEG 


Variable 

Description 

IT, JT 

Column and row indices for cell under consideration. 

SXX 


SYY 

STT 

Deviatoric components of stress. 

TXY 



Element Description: VISCOS calculates artificial viscosity. 


Abort Conditions: None 

Messages to Printer: None 


122 








Element Name: VOLCEL 


Element Type: FORTRAN Real 

Function Subroutine 

Entry Points: VOLCEL 

Called By: PRNTIJ, REGION, HYPOEL, PHASEA, STRESS, STRFLO , 

MASFLO, VARIBL 

External References: GETCON, ARECEL, SAVSUB, GETIJ, ABORT 

Common Blocks: CONSTS, SUBSBK, CUTBLK, RCDATA 

Calling Sequence: L, MODE 


Processor: 

I.C. , INTEG, EDIT 

Variable 

Description 

L 

Subscript to cell dictionary. 

MODE 

= 1; For regular cell . 
= 2; For lattice cell. 


Element Description: VOLCEL calculates the volume of a lattice or 

regular cell depending on mode . 

Abort Conditions: VOLCEL - The cell under consideration was cut but 

but there' was no location assigned to 
retrieve the area or volume from CUTBLK. 


Messages to Printer: None 
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Element Name: YLD 


Element Type: FORTRAN Real 

Function Subroutine 

Entry Points: YLD 



Called By: PIASFL 


External References: 

None 


Common Blocks: 

EOSBLK 


Calling Sequence: 

SF, EP 


Processor: 

INTEG 



Variable 

Description 

SF 

Logical variable which is true if work hardening is 


used. 

EP 

Plastic strain. 


Element Description: YLD provides the yield criterion when the von Mises 

criterion is not being used. 

Abort Conditions: None 

Messages to Printer: None 
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3.3 COLLECTOR SOURCE ELEMENTS . . . 

There are three COLLECTOR source elements in STEEP32. These are 
MAIN1, MAIN2; and MAIN3. These elements determine the memory alloca- 
tion structure for each of the three main processors of STEEP32. The 
COLLECTOR source language is described in Section 6.6 of the Marshall 
Space Flight Center Programmer's Procedure Manual. The primary output 
from the COLLECTOR is the absolute or relocatable element which will be 
executed on the UNIVAC 1108. It should be noted that the structure defined 
by the COLLECTOR also determines the amount of core required by a program. 
With judicious modification of the memory allocation scheme for a particular 
STEEP32 problem it might in some cases be possible to reduce the amount of 
storage required. 

The three diagrams on the following pages indicate the memory 
allocation scheme for the three main processors of STEEP32. 
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MEMORY ALLOCATION SCHEME FOR THE INTEGRATION PROCESSOR 



126 










MEMORY ALLOCATION SCHEME FOR THE INTEGRATION PROCESSOR (Continued) 
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MEMORY ALLOCATION SCHEME FOR THE INITIAL CONDITIONS PROCESSOR 
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MEMORY ALLOCATION SCHEME FOR EDIT PROCESSOR 
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SECTION 4 


STEEP 3 2 COMMON BLOCKS 


In this section each common block is described. A list of subroutines 
referencing the common block is given to aid any user desiring to make code 
changes. Finally, an extensive description of the variables in each common 
block appears. Several very important conventions involving the cell 
dictionary, the association of cut cell information and the location of the 
cell variables are described here. 
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COMMON BLOCK BOUNDY 


Used By: BDINIT, BDNORM , BDPOST, BDSTRS, BONDPR, COLDTS, CONTRL 

DATA, FILE, GENBND , INTPSH, MOVER, PLASFL, PLTLNT, 

PNTSET, PRTBPT, PUTGET, STRESS. 


Variable 

ENDSTR 


STRTSW 


N BOUND 


Description 

A logical variable which is set to true 
if the last point of a string is being 
processed. 

A logical variable which is set to true 
if the first point of a string is being 
processed. 

The total number of cut cells. NBOUND 
is used by FILE and INTPSH. 


NEXT Is the next available location in the 

push down list. NEXT is set by PUTGET 


POS Logical variable which is true when a 

real string is being processed. 

NSTRG The number of string in the problem. 

MSTRG (NUMSTG) The number of string points in each 

string . 

XSTRG (NUMSTP) The x-coordinates of the string points. 


YSTRG (NUMSTP) 


The y-coordinates of the string points. 


LSTRG (NUMSTP) 


The dictionary or control word for each 
string point. 
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COMMON BLOCK C BLOCK 


Used By: COLD, CONTRL, DUMSUB, DMSUB1, FIDE, GETCON, MASFLO , 

PHASEA, PRNTCT, PRNTIJ, PUTCON, SETDEL, STRESS, STRFLO 

C BLOCK contains only one array, the cell dictionary array - C(NUMCEL). 
The cell dictionary has one word for each cell of a STEEP32 problem. Into 
this word certain information about the cell is packed. TO extract 
information from the cell dictionary it is necessary to mask the word and 
shift the result. For this purpose the constant arrays MASK and SHIFT are 
used. The subscripts of MASK and SHIFT have been given names according 
to the following convention: 



C = regular cell 
L = lattice cell 
S = segment 
H = horizontal 
V = vertical 


Variable 

Name 

Number of 
Bits 

Description 

CMAT 

2 

Cell material 0 = empty, 1 = part 
filled, 2 = filled. 

LMAT 

2 

Lattice material 0 = empty, 1 = part 
filled, 2 = filled. 

CHS* 

3 

Cell horizontal segment type. 

CVS* 

3 

Cell vertical segment type. 

LHS* 

3 

Lattice horizontal segment type. 

LVS* 

3 

Lattice vertical segment type. 

LPT 

1 

Lattice point at n + 1 , 0 = outside 
of mass, 1 = inside of mass. 

CN . 

1 

Cell atn, 0 = no mass, l = mass. 

CNP 

1 

Cell at n + 1, 0 = no mass, 1 = mass 

LN 

1 

Lattice point at n, 0 = outside, 
1 = inside . 

LNP 

1 

Not used. 
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Common Block C BLOCK (concluded) . 

Variable Number of Description 

Name Bits 

LOC 10 Data location. If cell is cut this 

field will contain the location in 
CUTBLK to use to obtain the area 
and volume information. 


♦Segment types serve to establish which part of the cell segments contain 
mass. The meanings, of the different values are as follows: 


Value 

Definition 

0 

No mass 

1 

Cell mass 

2 

Mass, no mass 

3 

No mass , mass 

4 

No mass, mass, no mass 

5 

Mass, no mass, mass 
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COMMON BLOCK CON STS 


\ 


Used By: ARECEL, BDINIT, BDMOVE, BDPOST, BONDPR , CBN DRY, COLD, 

COMP1, COMP2, COMPUT, CONTRL, CONTUR, DATA, DATAL, 
DIVU, DMSUB1, DUMSUB, EXIST, EXTRAP, FILE, FILLER, 
FORDIC, GENBND, GETCON, GRIDVT , HYPOEL, LNEDIT, 

MAIN1 , MAIN 2 , MAIN 3 , MAPMIX, MASFLO, MASSIJ, PAPER , 
PAPER 1 , PHASEA, PIASFL, PLTSTR, PRNTCT , PRNTIJ, PUTCON, 
REGION, RESULT, SEGMNT, SELDT, SETDEL, SETFOR, SETMAT, 
SETSEG, SETSUB, SIO, STARVL , STRESS, STRFLO, TITLE, 

VARIBL, VECTOR, VISCOS, VOLCEL 


CONSTS contains most of the undimensioned constant variables. The majority 
of these variables are set by the namelist INITAL and are defined in the 
description of that namelist. The exceptions fall into two categories. CMAT, 
LMAT, CHS, CVS, LHS , LVS, LPT, CN, CNP, LN, LNP, LOC , MASK, and 
SHIFT all are used in conjunction with the cell dictionary word and are 
described with C BLOCK. STOPOP, STOPSM, STOPDT, STOPCY, STOPTI, 
STOPRT, STOPFO , STOPIL, STOPIR, STOPJL, STOPJU indicate the conditions 
which caused the problem to stop. These logical variables are set to true 
under the following conditions; 


Variables 

Description 

STOPOP 

The operator has set sense switch 5 
down. This is no longer used. 

STOPSM 

The value of the variable SMAX exceeds 
the maximum value allowed. 

STOPDT 

DT is less than DTMIN. 

STOPCY 

Current cycle exceeds MAXCYC-1. 

STOPTI 

Current time exceeds TMAX. 

STOPRT 

Not used . 

STOPFO 

The string has folded. 

STOPIL 

Motion on right hand boundary was 
sensed and SPRITE was true. 

STOPIR 

Motion on left hand boundary was 
sensed and SPLEFT was true. 

STOPJL 

Motion on bottom boundary was 
sensed and SPLOWR was true. 

STOPJU 

Motion on upper boundary was 
sensed and SPUPPR was true. 
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COMMON BLOCK CUTBLK 


Used By: ARECEL, BDPOST, COLD, CONTRL, FILE, GENBND, HYPOEL, 

MASFLO, PR^SEA, PIASFL, PRNTIJ, PRTBPT , PUTGET, SEGMNT, 
SETDEL, SETFOR, SETMAT, SETSEG, SETUP, STARVL , STPOST, 
STRFLO, VOLCEL. 


CUTBLK contains all the data for cut cells. Whenever a cell is cut, special 
information must be retained. Part of the information appears in the 
dictionary as segment types. The remaining data is referenced using the field 
LOC in the dictionary. If the dictionary is nonzero, either the regular cell or 
the stress cell is cut and the following information is recorded in CUTBLK: 
AREAC , AREAL, VOLS, VOLL, AL, RA, RB . 
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AREAC \ . f Area \ , /R , . % 

VOLC J 1S ^ Volume^ of ^1 + A 2 ) 

AREAL \ . /Area \ ' A 

VOLL f Volume j of A 3 


RA 


P^ is the first intersection on vertical segment. 


RB 


P 2 is second intersection on vertical segment. 


RA; 

RA 
AL 

for horizontal segment. At P]^ and ?2 since the 
vertical segment is cut two time the outer normal is 
not used, although one may be stored (depending on 
the ordering of the calculations the angle may be that 
associated with Pj or P 2 ). 


Pg is only intersection oh horizontal segment. 

S 1 is only intersection on horizontal segment. 

At Po the outer normal anale to the surface is stored 


Points P and are associated with adjacent cells . 
00 


136 


COMMON BLOCK DELSBK 


Used By: COLD, CONTRL, CONTUR, HYPOEL, MASFLO, PHASEA, PRNTIJ, 

SETDEL, SETFOR, STRESS, STRFLO 


DELSBK contains the communication data used with the flow of deviatoric stress 


Variable 


Description 


SXXH 

SYYH 

STTH 

TXYH 


Flow terms across the horizontal face. 


SXXV 

SYYV 

STTV 

TXYV 


Flow terms across the vertical surface. 


NOBDPT 


The number of cuts in cell sides. 


SINT 

COST 


Not used. 
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COMMON BLOCK DIVUBK 

Used By: DlVU, HYPOEL, PLASPL, STRESS, VISCOS 

DIVUBK contains the communication data associated with the calculation of 
the divergence. 


Variable 

DXX 

DXY 

DYY 

DTT 

DX 

DY 

Omega 

DMM 

partvl 



d 


zz 


d ee 

Length Of cell side in x direction. 
Length of cell side in y direction. 
Spin 
d 

mm 


Particle velocity. 
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COMMON BLOCK D STARS 
Used By: STRESS , VISCOS 


DSTARS contains the communication information between VISCQS and STRESS, 


Variable 


DTTS 


Description 
^00 ^00 ^mm 


DXXS 

DXYS 

DYYS 


d = d - d 
rr rr mm 



d* = d - d 
zz zz mm 
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COMMON BLOCK EOSBLK 


Used By: PRESS, REGION, STRESS, TBFT02 , VISCOS, YLD 


EOSBLK contains the communications data for the equations of state 
computations . 


Varia ble 
K 

PRES 

VL 

ETAL 

EL 


Description 

Pointer to the first PC ON constant. 
Pressure 

The sound speed squared. 

Relative density at lattice point. 

Specific internal energy at the lattice 
point. 
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COMMON BLOCK FRSTBK 


Used By: COMPUT, MAIN 2 


FRSTBK contains the logical variable FIRST. FIRST is set to true when the 
first cycle following a refresh is being integrated. This information is used 
to suppress the calculation of At for that. 
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COMMON BLOCK GRDATA 


Used By: BDMOVE, BOUNST, COLD, CONTRL, DMSUB1, DUMSUB, FILE, 

HYPOEL, LNEDIT , MASFLO , PHASEA, PLTSTR, PRNTIJ, REGION, 

• RESULT, SETDEL, SETFOR, SIO, STRESS, STRFLO, VARIBL, 
VECTOR, VISCOS 


GRDATA contains the core buffers for the cell variables. The names of the 
variables which are contained in GDATA depend on the stage of the calculation 
one is looking at. The following table establishes the relationship between 
the various cell variables and the subscripts of GRDATA. 


GRDATA 

Pass 1 

Pass 2 

Location 

COMP 1 Read 
(Values at n) 

COMP1 

COMP2 

Write 

Read 

COMP 2 Write 
(Values at n + 1) 


Value 

Value 

By 

Va lue 

By 

1 

E 

ETOT 

PHASEA 

E 

MASFLO 

2 

MASS 

MASS 

— 

MASS 

MASFLO 

3 

TTT 

STTA 

HYPOEL 

TTT 

STRESS 

4 

TXX 

SXXA 

HYPOEL 

TXX 

STRESS 

5 

TXY 

TXYA 

HYPOEL 

TXY 

STRESS 

6 

TYY 

SYYA 

HYPOEL 

TYY 

STRESS 

7 

U 

UAVE 

PHASEA 

U 

MASFLO 

8 

V 

VAVE 

PHASEA 

V 

MASFLO 

9 

DSXX 

STTEST 

HYPOEL 

DSXX 

STRESS 

10 

DSYY 

SXXEST 

HYPOEL 

DSYY 

STRESS 

11 

DSTT 

TXYEST 

HYPOEL 

DSTT 

STRESS 

12 

DSXY 

SYYEST 

HYPOEL 

DSXY 

STRESS 

13 

— 

UEST 

PHASEA 

— 

— 

14 



VEST 

PHASEA 

— 

— 

15 

— 

EEST 

PHASEA 

— 

— 

16 


TH 

PHASEA 
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COMMON BLOCK IOCNTL 


Used By: ABORT , BDMOVE , BONDPR, COLPTS, CONTRL, DATA, DATAL, 

FILE, GENBND , GRAPHS, HYPOEL, LEDPLT, LNEDIT , MAIN3, 
MAPMIX, MASFLO, PHASEA, PLOTTR, PRNTCT, PRNTIJ, PRTBPT, 
REGION, RESTIN, RESTOU, RESULT, SELDT, SEGMNT, SETDEL, 
SETFOR, SIO, STRESS, STRFLO , TITLE 


IOCNTL contains all logical unit assignments. 


Variable 

INP 

OUT 

BIN 

BOU 

PLT 

SIN 

SOU 

FRESH 


Description 

Input unit corresponding to card reader. 
Currently has value 5. 

Output unit corresponding to line 
printer. Currently has value 6. 

Unit for the input restart tape. 

Unit for the output restart tape. 

Unit for the SC4020 plot tape if any. 

Scratch drum file for read . 

Scratch drum file for write. 

Logical variable which is set to true 
when machine has been refreshed. 
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COMMON BLOCK MATBLK 

Used By: COM PUT, CONTRL, DATA, DATAL, FILE, HYPOEL, LOADRG , 

PIASFL, PRESS, REGION, STRESS, STRFLO 


MATBLK contains the material defining data . 


Variable 

Description 

TWOMU 

2 • ^ where p. is the shear modulus . 

KTYPE 

Not used. 

YIELD 

Coefficients to use in yielding. 

LFIT 

Equation of state number. 

LPOINT 

Pointer to location in PCON table 
to start. 

VONMIS 

A logical variable which is set to true 
by the user if he desires the Von Mises 
yield criterion to be employed. 
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COMMON BLOCK NETBLK 


Used By: CONTRL, DATA, DMSUB1, DUMSUB, GRIDVT, LOADRG 


NETBLK contains all the information necessary to generate the grid. 


Variable 

NOX 

NOY 

NEWXYS 

X(I) 

Y(J) 

DX(I) 

DY(J) 


Description 

The number of x's used to specify grid 

The number of y's used to specify grid 

Logical variable not currently used. 

x values to start equal spacing. 

y values to start equal spacing. 

Value to use to generate grid between 
X(I) and X(I + 1). 

Value to use to generate grid between 
Y(J) and Y(J + 1). 
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COMMON BLOCK PCNBLK 


Used By: DATA, DATAL, DMSUB1, DUMSUB, GRIDVT, LOADRG 


PCNBLK contains all the equation of state parameters. 


Variable 


Description 


PCON Array containing the equation of 

state parameters. The interpretation 
of the subscripts for PCON vary 
depending on the equation of state 
used . 
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COMMON BLOCK RCDATA 


Used By: ARECEL, BDMOVE , BDNORM , BDPOST, BDSTRS, COLD, COLPTS, 

CONTRL, CONTUR, DATA, DIVU, FILE, GENBND, GRIDVT, 
HYPOEL, LNEDIT , MASFLO , MASSIJ, PHASEA, PLASFL, PLTSTR, 
PRNTIJ , REGION, RESULT, SEGMNT, SETDEL, SETFOR, STRESS, 
STRFLO, VECTOR, VISCOS, VOLCEL 


RCDATA contains all the row and column vectors describing the problem grid. 


Varia ble 

Equivalent Name 

Description 

RDATA(1 , 1) 

YCS(l) 

y coordinate of regular 
cell side. 

RDATA(1 , 2) 

YCM(l) 

y coordinate of stress 
cell side. 

RDATA(1 , 3) 

DYCS(l) 

Length of regular cell side 
in y direction. 

RDATA(1 ,4) 

DYCM(l) 

Length of stress cell side 
in y direction. 

RDATA(1 , 5) 

ARCS(l) 

Cross-sectional area for 
regular cell . 

RDATA(1 , 6) 

ARCM(l) 

t 

Cross-sectional area for 
stress cell. 

CDATA(1 , 1) 

XCS(l) 

x coordinate of regular 
cell side. 

CDATA(1 , 2) 

XCM(l) 

x coordinate of stress 
cell side . 

CDATA(1 , 3) 

DXCS(l) 

Length of regular cell side 
in x direction. 

CDATA(1 , 4) 

DXCM(l) 

Length of stress cell side 
in x direction. 
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Association of cells, rows, columns, and segments (core file) RCDAIA 


XCM(I) XCM(I + 1) 



Subscription columns (I, LI) and rows (J, LJ) 


I takes on values 

IMIN 

thru 

IMAX 

J takes on values 

JMIN 

thru 

JMAX 

LI takes on values 

LIMIN 

thru 

UMAX 

LJ takes on values 

LJMIN 

thru 

LJMAX 


Given I and J, or LI and LJ, subroutine SETSUB will compute IJ and/or LIJ 
and IJD. 

IJ and LIJ are addresses in the grid buffers (GRDATA) . IJD is adress in 
dictionary . 
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COMMON BLOCK RCTEMP 


Used By: CONTRL, CONTUR, MAPMIX, MASFLO, PHASEA, PRNTIJ, SETDEL, 

SETFOR, STRFLO 


RCTEMP contains the row and column scratch information. 


Variable 

Equivalent Name 

Description 

RTEMP(1 , 1) 

DEH(l) 

The change in energy across the 
horizontal surface of a cell. 

RTEMP(1 , 2) 

DEV(l) 

The change in energy across the 
vertical surface of a cell. 

RTEMP(1 , 3) 

DMH(l) 

The change in mass across the 
horizontal surface of a cell. 

RTEMP(1 ,4) 

DMV(l) 

The change in mass across the 
vertical surface of a cell. 

RTEMP(1 , 5) 

DUH(l) 

The change in U momentum 
across the horizontal surface 
of a cell. 

RTEMP(1 , 6) 

DUV(l) 

The change in U momentum 
across the vertical surface of 
a cell. 

RTEMP(1 , 7) 

DVH(l) 

The change in V momentum 
across the horizontal surface 
of a cell . 

RTEMP(1 , 8) 

DW(1) 

The change in V momentum 
across the vertical surface of 
a cell. 



" COMMON BLOCK REGBLK ' 

Used By: CONTRL, DATA, DMSUB1, DUMSUB, LOADRG , REGION 


REGBLK contains the region specifications which are used to define the initial 
conditions. The definition of the variables used in REGBLK can be found with 
the namelist INITAL . The variables in REGBLK are': NREG, XLREG , XRREG , 
YLREG, YUREG, UREG, VREG, EREG, RHOREG, STTREG, SXXREG, SYYREG, 
TXYREG, REGTYP. . ... . 


COMMON BLOCK SCANBK 
Used By: BDPOST, CONTRL, SETSID 


SCANBK contains the communication data for the dictionary manipulation. 


Variable 

MD 

MA 

MB 

ME 

LA 

LB 

LE 

KEY 


Description 

The dictionary subscript for the cell. 

The minimum column number minus 1 . 

The maximum column number. 

The column offset for the dictionary. 

The minimum row number. 

The maximum row number. 

The row offset for the dictionary. 

The key to the segment being 
considered. 
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COMMON BLOCK SCCNTL 

Used By: CONTRL, CONTUR, FUNCT, GRAPHS, LEDPLT, LNEDIT , PAPER, 

PAPER 1 , PLOTTR, PLTBND , PLTSTR, RESULT, VECTOR 


SCCNTL contains all the control information for the SC4020 edits. The 
relevant variables are documented with the namelist SC4020. 


COMMON BLOCK SCDEPS 

Used By: LEDPLT, LNEDIT, PAPER, PAPER1, PLOTTR, RESULT 


SCDEPS contains all the variables connected with the linear editor. These 
variables are discussed in the linear editor documentation. 
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------ COMMON BLOCK SEGBLK ' 

Used By: BOUNST, COLD, CONTRL, CONTUR, HYPOEL, MASFLO , PHASEA, 

PRNTIJ, SEGMNT, SETDEL, SETFOR, STRESS, STRFLO 


SEGBLK contains information relating 

Variable 

NOSEG 

NOCUT 

XA 

XB 


to the segment computations. 
Description 

The number of segments in cell side 
under consideration. 

The number of cuts. 

y coordinate of intersection on a 
vertical segment. 

x coordinate of intersection on a 
horizontal segment. 
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COMMON BLOCK SHCOM 


Used By: COLD, COMPUT, CONTRL, DATA, FILE , GENBND, PLASFL, 

PRTBPT, STPOST, STRESS 


SHCOM contains the communication data used in work hardening and the 
Perzyna model. 


Variable Description 

SF Logical variable which is true if work 

hardening is to be performed. 

NIM The number of imaginary string points 

in problem . 

XBCOMP x coordinate of left side of column. 


XTCOMP x coordinate of right side of column. 

DXCRAT Distances in x direction from the 

center of the current column to each 
of the imaginary string points. 

LPNT An array used by PLASFL to temporarily 

store subscripts of the imaginary 
string points which are being used for 
the work hardening calculations. 
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COMMON BLOCK SIOBLK 


Used By: MAIN1, MAIN 2 , SIO 

SIOBLK contains three logical variables. 


Variable 

Description 

CLOSED 

Logical variable which is set to true 
if the file has been closed. 

OPEN 

Logical variable which is set to true 
if the file has been opened. 

NEVER 

Logical variable which has the value 
true if the file has never been 


written. 
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SECTION 5 


INPUT DATA TO STEEP32 


Most input variables to STEEP32 occur in one or more namelist. r 
A particular namelist may be referenced by more than one processor. The 
sequence of the input data for each processor is described below. This 
is followed by a detailed description of each of the namelists. 


5.1 SEQUENCE OF INPUT DATA 

5.1.1 Initial Conditions 
$INITAL 
• - 

data for the INITAL namelist 


$END 

Heading card 
STARTR card 
$BDSTRG 


string data for problem 


$END 

CUTOFF card. 

The heading card is a card with an alphanumeric description of the 
problem in columns 1 to 72 . The STARTR card has the word STARTR in 
columns 11 to 16. The CUTOFF card has the word CUTOFF in columns 11 

to 16. 


5.1.2 Integration 
Tape Card 
$DATAL2 


specific data for problem 


$END 

Heading card if the logical variable NEWHED is set to true. 
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The format of the tape card is 514 . The variables appearing on it 
are BIN, BOT, SIN, SOU, and INFILE. These variables are the restart 
input tape, the restart output tape, the two scratch drum units, and the 
file number to use for restart. As with the initial conditions, the heading 
card contains 72 columns of alphanumeric information. 


5.1.3 Edit 
$INITAL 

m 

BIN, SIN, and SOU must be set here 

• 

$END 

EDIT card 
$SC4020 

data for plots 


NGRAPH pairs of these must appear if the 
logical variable NEWHED is true. 

The EDIT card must have the word EDIT appearing in columns 11 
through 14. 

$RSTART 

INFILE must be set here 
$END 


$END 

Heading Card 
Label Card 
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c 

r* 

5.2 NAMELIST DESCRIPTIONS 



U 

c 

c 

5.2.1 THE 

NAMELIST INITAL 


\ 

c 

c 

the variables for 

the NAMELIST INITAL ARE DIVIDED 

INTO GROUPS 

c 

according 

to their function. 



c 

within each group 

THE DESCRIPTION OCCURS IN FOUR ' 

COLUMNS. THE 

c 

first column gives the variable name, the second column gives 

c 

the variable description, the third column .indicates the variable 

c 

TYPE, the 

FOURTH 

COLUMN DEFINES THE ASSUMED VALUE 

THE VARIABLE IS 

c 

assigned. 





c 

control data 





c 

BIN 


INPUT RESTART TAPE 

I 

0 

c 

BOT 


output restart tape 

I 

0 

c 

cycle 


cycle number 

I 

0 

c 

DT 


INTEGRATION step 

R 

. 1 

c 

DTMIN 


MINIMUM OT 

R 

.01 

c 

FRESH 


MACHINE STATUS. 

L 

.false. 

c 



if fresh is true, then there 



c 



is a good file on drum file 



c 



sin. 



c 

MAXCYC 


MAXIMUM NUMBER INT. CYCLES 

I 

200 

c 

nbound 


MAXIMUM NO. BOUNDARY CELLS 

I 

100 

c 

NP A USE 


OUTPUT FREQUENCY 

I 

1 

c 

SF 


true for strain hardening 

L 

.false. 

c 

SIN 


scratch input file 

I 

15 

c 

SOU 


scratch output file 

I 

16 

c 

TIME 


time 

R 

0. 

c 

TMAX 


maximum int. TIME 

R 

0. 

c 

NEWHEO 


HEADER card read control 

L 

.false. 

c 

GEOMETRY AND 

boundary 

DATA 



c 

GEOM 


GEOMETRY 1 =sl ab » ?=CYL. 

I 

0 

c 

LBRITE 


RITE BOUNDARY l=CONT. ,2=FIXED 

I 

i 

c 

LBLOWR 


LOWR BOUNDARY l=CONT. .2=FIXED 

I 

i 

c 

LRLEFT 


LEFT BOUNDARY 1 =CONT . » 2=F I XED 

1 

l 

c 

L6UPPR 


UPPR BOUNDARY I =CONT . . 2=F I XED 

I 

l 

c 

SPLEFT 


STOP FOR MOTION LEFT BNDY. 

L 

.FALSE. 

c 

sprite 


STOP FOR MOTION RIGHT BNQY. 

L 

.false. 

c 

SPLOWR 


STOP FOR MOTION LOWER BNDY. 

L 

.false. 

c 

SPUPPR 


STOP FOR MOTION UPPER BNDY. 

L 

.false. 

c 

NCOL 


NUMBER OF CFLL COLUMNS 

I 

i 

c 

NROW 


NUMBER OF CELL ROWS 

I 

l 

c 

NOX 


NUMBER OF X.DX COMBOS 

I 

0 

c 

X(l) 


LEFT X«S 

R 

0. 

c 

DX ( 1 ) 


DX INCREMENTS 

R 

0. 

c 

NOY 


NUMBER OF Y.DY COMBOS 

I 

0 

c 

Y( 1) 


LOWER Y*S 

R 

0. 

c 

DY ( 1) 


DY INCREMENTS 

R 

0. 

c 

NREG 


NUMBER OF REGIONS 

I 

0 

c 

XLREG(l) 


LEFTMOST X OF REGION 

R 

0. 

c 

XRREG (1) 


RIGHTMOST X OF REGION 

R 

0 . 

c 

YLREG ( 1 ) 


LOWERMOST Y OF REGION 

R 

0 . 

c 

YUREG ( 1 ) 


uppermost y of region 

R 

0 . 

c 

EREG ( 1 ) 


specific internal volumetric 

R 

0 . 

c 



ENERGY OF REGION 



c 

UREG ( 1 ) 


x VELOCITY OF REGION 

R 

0 . 

c 

VREG(l) 


Y VELOCITY OF REGION 

R 

0 . 

c 

RHOREG ( I ) 



R 

0 . 

c 

TTTREG(l) 


TTT STRESS OF REGION 

R 

0 . 
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c 

TXXREGtl ) 

TXX STRFSS OF REGION 


R 

c 

TXYREG ( 1 ) 

T X Y STRFSS OF REGION 


R 

c 

TYYREG ( 1 ) 

TYY STRFSS OF REGION 


R 

c 

REGTYP(I) 

TYPE OF REGION ( 0 = ST ANDARD) 

I 

c 1 

MATERIAL DATA 




c 

TWOMU 

Z* MU' FOP STRENGTH EFFECTS 

R 

c 

VONMIS 

VON meisis criterion 


L 

c 

YIELD 

YIELD CONSTRAINT. IF=0. ♦ 

NONE 

R 

C 1 

EQUATION OF STATE 

DATA 



c 

LFIT(I) 

EOS FORM . 


I 

c 

LPOINT (1) 

POINTER TO EOS DATA 


I 

c 

PCOM(l) 

EOS DATA 


R 

c 

PCONTe 

CONTROL FOR READING PCON 

data 

L 

c 

STABILITY data 




c 

Si 

stability control si 


R 

c 

S2 

STABILITY CONTROL S2 


R 

c 

R 1 

INCREASE FACTOR FOR DT 


R 

c 

R2 

DECREASF FACTOR FOR DT 


R 

c 

TOi FRANCES 




c 

PTOL 

MINIMUM ABSOLUTt PRESSURE 

R 

c 

DUTOL 

TOLERANCE ON QU (aBS) 


R 

c 

DVTOL 

TOLERANCE ON f)V (ABS) 


R 

c 

OXTOL 

MINIMUM OX FOR BOUNDARY 

POINT 

R 

c 

DYTOL 

MINIMUM OY FOR BOUNDARY 

POINT 

R 

c 

ETOL 

TOLERANCE ON E (ABS) 


R 

c 

RHOTOL 

TOLERANCE ON RHO (ABS) 


R 

c 

stttol 

TOLERANCE ON TTT (ABS) 


R 

c 

SXXTOL 

TOLERANCE ON TXX (ABS) 


R 

c 

SXYTOL 

TOLERANCE ON TXY (ABS) 


R 

c 

SYYTOL 

TOLERANCE ON TYY (ABS) 


R 

c 

UTOL 

TOLERANCE ON U (ABS) 


R 

c 

voltol 

TOLERANCE on rel. volume 

R 

c 


FOR ACCELERATION, 

WORK 


c 

VTOL 

TOLERANCE ON V (ABS) 


R 


C VISCOSITY DATA 

C COEVIS(l) linear* quadratic, DEVIATOR 

c VISCOSITY COEFFICIENTS 

c QMAXU) (UNUSED) 

c$ 
c 


0 . 

0 . 

0 . 

0 

0 

.TRUE. 

0 . 

0 

1 

0. 

F 

.16 

.35 

1.5 

.5 

l.E-6 
l.E-6 
l.E-6 
.001 
.001 
l.E-6 
l.E-6 
l.E-6 
l.E-6 
l.E-6 
1 .E-6 
l.E-6 
.1 

l.E-6 

0. 


158 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

( 51(3 


5.2.2 The namelist bdstrg 

THE NAMELIST 8DSTRG IS USED TO DEFINE THE STRING POINTS. 


VARIABLE _ DESCRIPTION 

name 


TYPE ... ASSUMED. 
VALUE 


EPS 

NIM 


nstrg 

MSTRG ( i ) 


XSTRG ( 1 ) 
YSTRG(l) 


INIT.AL EPS STRAIN HARDENING R o. 

IE (SF)NIJMBER OF STRAIN POINTS I 0 

STRAIN POINTS MUST PRECEDE 
EITHER OTHER IMaGINAY POINTS 
OR OTHER REAL STRING PNTS. 

NIM MUST PE .GE. A 

NUMBER OF STRINGS I 0 

NUMRER OF POINTS PER STRING I 0 

IF NEGATIVE* DUMMY STRING 
OR STRAIN hardening POINTS 
X COORDINATE of STRING POINT R 0. 

Y COORDINATE OF STRING POINT R 0. 


NOTE THAT THERE ARE TWO TYPES OF STRINGS* IMAGINARY AND REAL. THE 
REAL STRINGS DEFINE MATERIAL. AS THE SUBCRIPTS OF THE STRING 

points increase* the material lies to the right, since steep32 is 
a single material code* this implies that there is void on the 
left of the string. IMAGINARY string points do NOT DEFINE MATERIAL 
AND AS SUCH they HaVE NO FEFECT On THE AREA AMD VOLUME 
CALCULATIONS, imaginary string serve FIRST OF all to trace 
MATERIAL FLOW, in PROBLEM USING STRAIN HARDENING THEY MAY 
ALSO 8E USED TO CARRY THE PLASTIC STRAINS. 


5.2.3 THE NAMELIST rstapt 
THE NAMELIST RSTaRT IS USED TO 

the eoit processor is plotting 


CONTROL the RESTART FILE THAT 
FROM. n 


VARIABLE DESCRIPTION TYPE ASSUMED 

NAME VALUE 

INFILF FILE NUMBER FOR REFRESHING I 0. 


*1 
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5.2.4 THE NAMELIST DAT AL2 


C 

c 

c The NAMELIST DaTaL2 is read by the integration processor to 

c CHANGE VARIABLES FROM THOSE VALUES THEY HAD ON THE RESTART TAPE. 

C SOME OF THE VARIABLES SUCH AS MaXCYC WILL GENERALLY REQUIRE 

C RESETTING WITH EACH RESTART. 

C 

c 


c 

name 



VALUE 

c 

variable 

DESCRIPTION type 

ASSUMED 

c 

c 

DT 

INTEGRATION step 

R 

.1 

c 

NPAUSE 

OUTPUT FREQUENCY 

I 

1 

c 

MAXCYC 

MAXIMUM NUMBER INT. CYCLES 

I 

200 

c 

tmax 

MAXIMUM INT. TIME 

R 

0. 

c 

DTMIN 

MINIMUM OT 

R 

.01 

c 

SI 

stability control si 

R 

.16 

c 

S2 

stability control s? 

R 

.35 

c 

R 1 

increase factor for dt 

R 

1.5 

c 

R2 

decrease factor for dt 

R 

.5 

c 

DUTOL 

TOLERANCE ON DU (ABS) 

R 

l.E-6 

c 

DVTOL 

TOLERANCE ON QV (ABS) 

R 

1 .E-6 

c 

TXXTOL 

TOLERANCE on TXX (ABS) 

R 

l.E-6 

c 

TXYTOL 

TOLERANCE ON TXY (ABS) 

R 

1 .E-6 

c 

TTTTOL 

TOLERANCE ON TTT (ABS) 

R 

1 .E-6 

c 

TYYTOL 

TOLERANCE ON TYY (ABS) 

R 

1 ..E-6 

c 

UTOL 

TOLERANCE ON U (ABS) 

R 

1 .E~6 

c 

VTOL 

TOLERANCE ON V (ABS) 

R 

l.E-6 

c 

ETOL 

TOLERANCE ON E (ABS) 

R 

l.E-6 

c 

RHOTGL 

TOLERANCE on RHO (ABS) 

R 

1 .E-6 

c 

TWOMU 

2*MU FOR STRENGTH EFFECTS 

R 

0. 

c 

LF IT ( 1 ) 

EOS FORM 

I 

0 

c 

LPOINT (1) 

POINTER TO EOS DATA 

I 

1 

c 

PCON ( 1 ) 

EOS DATA 

R 

0. 

c 

lbrite 

RITE BOUNDARY 1=C0NT. .2=FIXED 

I 

1 

c 

. lbleft 

LEFT BOUNDARY 1=C0NT . «2=FIXED 

I 

1 

c 

lblowr 

LOWR BOUNDARY 1=C0NT. *2=FIXED 

I 

1 

c 

LeUPPR 

UPPR BOUNDARY 1 =CONT . • 2=F I XtO 

I 

1 

c 

SPLEFT 

STOP FOR MOTION LEFT BNDY . 

L 

.false. 

c 

SPRITE 

STOP FOR MOTION RIGHT BNDY. 

L 

.false. 

c 

SPLOWR 

STOP FOR MOTION LOWER BNDY. 

L 

.false. 

c 

SPUPPR 

STOP FOR MOTION UPPER. BNDY. 

L 

.FALSE. 

c 

BIN 

INPUT RESTART TAPE 

I 

0 

c 

.80 T 

OUTPUT RESTART TAPE 

I 

0 

c 

SIN 

scratch input file 

I 

15 

c 

sou 

SCRATCH OUTPUT FILE 

I 

16 

c 

c 

voltol 

tolerance on pel. volume 

FOR ACCELERATION* WORK 

R 

. 1 

c 

0XT0L 

MINIMUM DX FOR BOUNDARY POINT 

R 

.001 

c 

0YT0L 

MINIMUM DY FOR BOUNDARY POINT 

R 

.001 

c 

NEWHED 

HEADER CARD READ CONTROL 

L 

.FALSE. 

c 

c 

COEVlS(l) 

linear* quadratic, oeviator 

VISCOSITY COEFFICIENTS 

R 

0. 

c 

PT0L 

MINIMUM ABSOLUTE PRESSURE 

R 

l.E-6 

c 

pcontb 

CONTROL FOR READING PCON DATA 

L 

F 

c 

V0NMIS 

VON MISES CRITERION 

L 

.TRUE. 

c 

(5)i3 

SF 

TRUE FOR STRAIN HARDENING 

L 

.FALSE. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


5.2.5 The NAMELIST SCA020 


THE NAMtllST SC4020 CONTAINS ALL The VARIABLES USED TO DEFINE AN 

SC4020 plot. 


variable 

NAME 

PAPERT(l) 


XMIN ( 1 ) 
XMAX(l) 

YM I N ( 1 ) 

YM AX ( 1 ) 
DELX (1 ) 
DELY ( 1 ) 
ALPHA 1 (1) 
ALPHA? ( 1 ) 
FCLASS (I) 


■» 

* 


XSCALE ( 1 ) 
TYPLEN( 1 ) 
YSCALE ( I ) 
FMIN ( 1 ) 
FSCALE ( 1 ) 
FUNCX ( 1 ) 
FUNCY ( 1 ) 

WHEN LNEDIT 
FUNCX ( 1 ) 

FUNCY (1) 


IS 


FUNCY ( 1 ) 


PLOT 
LINEOT 
NPEGON 
MINI(l) 
MAXI ( 1 ) 
MINJ ( l ) 
MAX J ( 1 ) 


description 


type assumed 

VALUE 
I 0 


paper type 

=1 EQUAL X,Y SCALE 
= ? standard GRID 
= ? ISOMETRIC 

MINIMUM X FOR GRAPH ' R 0. 

MAXIMUM X FOR GRAPH R . 0. 

MINIMUM Y FOR GRAPH . R 0. 

MAXIMUM Y FOR GRAPH R 0. 

DELTA X FOR GRAPH LINES R 0. 

DELTA Y FOR GRAPH LIMES R . 0. 

ALPHA/ 1 FOR ISOMETRIC R 0. 

ALPHA/? for isometric R 0. 

FUNCTION CLASSIFICATION I 0 

=1 FIX) VS F ( Y ) 

=2 velocity vectors 
=3 ISOMETRIC 
=4 MASS POSITIONS. 

=5 CONTOURS 

=6 PRINCIPAL STRESS 

=7 MAXIMUM SHEAR 

NOT IMPLEMENTED AS OF 03/20/72 


DIVISOR 
TYPICAL 
DIVISOR 
MINIMUM 
DIVISOR 
X FUNCTION 
Y FUNCTION 


FOR X FUNCTION R 0 

LENGTHS FOR PLOTS R 0 

FOR Y FUNCTION R 0 

length vector plotted r 0 

FOR VECTOR LENGTH R 0 

IN FX VS FY I 0 

IN FX VS FY I 0 


IMPLEMENTED FUNCX, FUNCY 

=i omega bar 
=? omega 

FOR FUNCX=] ONLY 
= 1 
= 7 
= 3 
= 4 
= 5 


take on following values 


MASS 

KINETIC ENERGY 
X -momentum 
y-momentum 

MOMENTUM MAGNITUDE 


CUMULATIVE 
CUMULATIVE 
CUMULAT I VE 
CUMULAT I VE 

cumulative 

FOR FUNCY=? ONLY 

= MASS PER STERADJAN 
=7 KINETIC ENERGY per steradian 
= o X-MQM PF-P STERADIAN 
=9 Y-MOM per STERADIAN 
=IOMOm mag per steradian 
PLOT COnTROL-IF F NO SCA020 L T 

linear FXTRaP editor control l f 

NO OF RFG SGMTS IN HONOR REG 1 0 

min j- value for donor reg i i o 

MAX I VALUE FOR DONOR REG I I 0 

MIN J Value FOR DONOR REG 1 I 0 

MAX J VALUE for DONOR REG 1 I 0 
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c 

XEDIT 

X VALUE FOR EXTRAPOLATION 

R 

0. 

c 

NSECTN 

NO OF SECTIONS+I AT XEDIT 

I 

0 

c 

SECTN(l) 

y-value of section divisions 

R 

0. 

c 

I TYPE 

lnedit PRINT CONTROL 

I 

0 

c 


=1 case l-XC.XMOM(J) 



c 


-2 CASE 2-XC i X-AVE 



c 


=3 CASE 3-XCtX-MASS-AVE 



c 


=4 all cases 



c 

CONI < I ) 

CONNECT EQUAL I IN FX VS FY 

I 

0 

c 

CONJ ( 1 ) 

CONNECT EQUAL J IN FX VS FY 

I 

0 

c 

NGRAPH 

NUMBER of graphs/file 

I 

0 

c 

nedit 

NUMBER of FILES TO EDIT 

I 

1 

c 

DOaLL 

CONTROL FOR ALL FILES/TAPE 

L 

F 

c 

CONTYP ( 1 ) 

type of contour to plot 

I 

0 

c 


=1 SPECIFIC INTERNAL energy 


c 


=? DENSITY 



c 


=3 DYNAMIC PRESSURE 



c 


=4 kinetic energy 



c 


=S MOMENTUM u-direction 



c 


=6 MOMENTUM v-direction 



c 

NCONT ( I ) 

NO. OF CONTOURS FOR EACH GRAPHI 

0 

c 

POINTR ( 1 ) 

pointer to ist value in table 

I 

I 

c 

SPACNG ( 1 ) 

TABLE OF CONTOUR VALUES 

R 

0. 

c 

PEW 

REWIND CONTROL riEFORE EDIT 

L 

T 

c 

IMINPT < 1 ) 

MIN I TO PLOT IN FX VS FY 

I 

0 

c 

IMAXPT ( 1 ) 

MAX I TO PLOT IN FX VS FY 

I 

0 

c 

JMINPT ( 1 ) 

MIN J TO PLOT IN FX VS FY 

I 

0 

c 

JMAXPT ( 1 ) 

MAX J TO PLOT IN FX VS FY 

I 

0 

c 

camera 

CAMERA 

I 

935 

c 

maxfrm 

MAXIMUM FRAME COUNT 

I 

50 

c 

numi 

(NOT USFD) 

I 

0 

c 

NUM? 

(NOT USED) 

1 

0 

c 

NEWHED 

CONTROL TO READ HEADING CARDS 

L 

T 

c 

LAST 

CONTROL FOR END OF PLOTS 

L 

F 

c 

ndecl 

THE MUMPER OF DECIMAL 

I 

2 

c 


PLACES TO THE LEFT OF THE 



c 


decimal point to be USED 



c 


IN LABELING THE X AND Y AXIS. 




0(5 
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SECTION 6 


A DESCRIPTION OF THE PRINTED OUTPUT FROM STEEP 3 2 


The printed output from STEEP32 falls into four categories. These 
are the printout of input data, cycle by cycle printout, restart file printout 
and diagnostic printout. 


6.1 INPUT DATA PRINTOUT 

The printout of input data generally takes the form of a WRITE 
statement immediately following the READ statement used to input the data 
Thus the format of this printout is the same as the input data described 
earlier . 


6.2 CYCLE BY CYCLE PRINTOUT 

An example of cycle by cycle printout appears on the next page. 
This group of print begins with the statement CYCLE=1 and continues 
through the next nine lines. The first line indicates the current cycle 
number and problem time. The next line contains totals for the mass, 
momenta, and energies in the problem. Lines three through six contain 
the peak stress and pressure as well as the indices of the cells where 
these peaks occur. 

In the example given, the cycle .1 printout is followed by two lines 
and the message "TOTAL NEGATIVE ENERGY SUPPRESSED." These three 
lines are produced by subroutine MASFLO and are not part of the regular 
printout. Following these three lines of diagnostic printout the cycle by 
cycle print for cycle 2 begins . 
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FILE WRITTEN ON EOT 2 6 NFILE 5 2 CYCLE 2 0 TIME- 0,0000 18 MAR 7* 10 * 07 1 45 OUTPUT 

COMPUT 


CYCLE" 1 TIME 5 1,2700-02 


MASS 2 

2.6998-01 umom= 


0,0000 


VMOM = 

5,7725-02 K.E, 5 


2, 

1264- 

02 I , E , 5 

3,8321-04 T , E , 2 

2 i 1647-02 

TXXMAX 2 

,4/600-01 AT 

LI 2 

5 LJ 2 

13 

TXXMIN 2 

- , 64QQ0*02 a T 

LI 5 

4 

LJ= 

12 



TYYMAX 2 

.16800-01 AT 

LI 2 

4 LJ 2 

12 

TYYMIN 2 

-.95200-01 AT 

LI 2 

5 

LJ = 

13 



TTTMAX 2 

,4/600-01 AT 

L I 2 

5 LJ = 

13 

TTTMIN 2 

-.84000-02 AT 

LI 2 

4 

LJ 2 

12 



PMAX = 

,ouooo at 

L I 2 

5 LU = 

37 

PM I N 2 

,00000 AT 

LI 2 

5 

LJ 5 

37 


i 


INCREASED DT TO 1 . 9050UQO-02 , SMAX 2,0250002-01 

DT= 1 , 9Q5U-U2 TIME 2 1,2700-02 CYCLE 2 1 Sl= 1,6000*01 S2s 3,50QO H 0l Rla 1,5000 + 00 R2? 5,0000«0i COMPUT PROCESSOR 

SMAX 2 2. .0250-01 CMaX 5 0,0000 VMAX 2 2,0250-01 ICMAX* 0 JCMAX= 0 IVMAX? 3 JVM'AX 5 12 CSMAX 2 O.OOoQ 

CYCLE 2 2 TIME 2 3,1750-02 

MASS 2 2 , 6VV8-01 UMOM= 0,0000 VMOM 2 5,7725-02 K.E, 2 2.0746-02 j.E, 5 9 ,qq51-o4 T,E,= 2i1647-o2: 

TOTAL NEGATIVE ENERGY SUPPRESSED I 


CYCLE 2 2 TIME 2 3,1750-02 


MASS 2 

1 , 7 158-03 UMCM 2 

0,0000 


V M 0 M 2 

5,8027-03 K 

>E,= 

2.1819 

"03 i,e,? 

-9,2047-08 T , E , 2 

2 1 1819-03 

TXXMAX 2 

,91792-01 AT LI 2 

3 LU = 

13 

TXXMIN? 

32301*01 

AT LI 2 

3, LJ 2 

12 > • 


i 

TYYMAX 2 

,64602-01 AT L I 2 

3 LJ = 

12 

TYYMIN 2 

-.18358-00 

AT LI 2 

5 LJ = 

13 


! 

TTTMAX 2 

,91792-01 AT LI 2 

3 U 2 

13 

TTTMIN 2 

-.32301-01 

AT LI 2 

3 LJ 5 

12 


1 

PMAX 2 

• OUOOO AT l I 2 

5 LJ 2 

37 

PMIN 2 

,00000 

AT LI 2 

5 LJ = 

37 


I 


DT= 1,905^02 TIME= J , 1750-02- CYCLE 5 -. 2 SI 2 1,6000-01 S2= 3,5000-01 Rl= 1.5000+00 R2= 5,0000*01 COMPUT PROCESSOR 

SMAX 2 2,0^15-01 CMAX- 0.0000 VMAX 2 2,0315-01 ICMAX 2 0 JCMAX 2 0 JVMAX 2 5 JVMAX 2 12 CSMAX 2 0,QboO 


i 

I 
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6.3 RESTART CYCLE PRINTOUT 

Whenever a restart file is written onto tape a detailed printout of the 
problem occurs. The three pages that follow are examples of restart file 
printout . 

The first page of printout is a display of the cell dictionary words. 
This is produced by subroutine PRNTCT. There are three groups of columns 
occurring on a page. Within a group the first two columns contain the cell 
indices. The dictionary subscripts are in the third column. Following this 
are twelve columns containing the values for the entries in the cell 
dictionary. For a description of the entries in the cell dictionary reference 
should be made to the documentation of common block CBLOCK in Section 4. 

The second page of output contains the locations of the strings and 
string points. The columns headed by EP and DELEP are empty since strain 
hardening was not being performed. 

All the grid and cell variable information appears on the third page 
of output. The first group of print contains the values of variables which 
are defined at the center of the regular cell. The second group of print 
contains the values of the lattice point variables. 

In addition, the restart cycle printout usually contains the printout 
produced by subroutine MASFLO. 


6.4 DIAGNOSTIC PRINTOUT 

Diagnostic printouts may occur from time to time during the integra- 
tion process. These diagnostics are intended to alert the user to possible 
abnormal situations. The messages produced by each subroutine are 
documented as part of the FORTRAN element description appearing in 
Section 3.2. 
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PURE ELAS, PRUB, USING OLD PCF FOR COMPARISON WITH SOLN FROM ST£EP32 SHOCK HYDRODYNAMICS,, INC, STRIPE C^DE 17 MAR 71 I3jlli01 
CYCL£= 0 T IM£= 0.0000 DT= 1.2700-02 STARTR PROCESSOR ...... 


I 

J 

I J 

mat 

CHV 

LHV 

L 

CNP 

LNP 

LOC 

4 

34 

0227 

2 

2 

1 

1 

1 

1 

1 

0 

1 

n 

U 

0 

000 

4 

5 7 

0232. 

2 

2 

1 

1 

1 

1 

1 

0 

1 

0 

0 

non 

5 

4 

X 

0235 

0 

0 

0 

n 

u 

0 

0 

0 

n 

0 

Q 

0 

non 

5 

4 

024 0 

1 

2 

1 

3 

1 

1 

1 

0 

1 

Q 

0 

454 

5 

7 

02 4 3 

2 

2 

1 

1 

1 

1 

1 

0 

1 

0 

0 

000 

5 

10 

0246 

2 

2 

1 

1 

1 

1 

1 

0 

1 

n 

0 

000 

5 

13 

0251 

2 

2 

1 

1 

1 

1 

1 

j 

1 

0 

0 

000 

5 

16 

0254 

2 

2 

«L 

1 

1 

1 


0 

1 

0 

a 

000 

5 

IV 

0257 

2 

2 

1 

4 

X 

1 

1 

i 

A 

0 

1 

0 

0 

oon 

5 

2 2 

0262 

2 

O 

1 

1 

1 

1 

1 

0 

1 

0 

0 

000 

w 

.y 

25 

0265 

2 

2 

1 

1 

4 

X 

1 

1 

0 

4 

X 

0 

0 

non 

6 

23 

0270 

2 

2 

1 

1 

1 

1 

1 

0 

1 

0 

0 

non 

5 

-* *, 

0273 

2 

2 

1 

1 

1 

1 

1 

0 

1 

0 

0 

.000 

5 

3 4 

027 6 

2 

2 

1 

1 

1 

1 

1 


1 

n 

0 

non 

5 

37 

C5Q1 

2 

2 

1 

1 

4 

X 

1 

1 

0 

1 

0 

0 

000 


10 IJ MAT CHV LHV L CNp ENP LOC 
4 35 0230 22111110100 000 

4 33 0233 j 1 020200 1.00440 

5 2 0236 0 0 0 0 0 0 0 0 C 0 0 000 

5 5 0741 22111110100 000 

5 8 0244 22111110100 000 

5 11 0247 22lllll0i 00 000 
5 14 Q232 221111101 OQ 000 
5 17 0235 2 2 1 1 1 1 1 0 1 o 0 000 
5 20 o 2 *o 2 2iiiiioinnoon 
5 23 02*3 2 2 1 1 1 1 1 0 l.o 0 000 
5 26 02*6 221111,10100 OOH 
5 29 0271 22111 1 1 0 1 q 0 000 
5 32 02 7 4 22111110100 000 
5 35 Q? 7 7 2 2 1 1 1 1 1 0 1 o 0 000 
5 38 0302 1 1 0 2 0 2 0 0 1 0 0 437 


I j u mat chv lhv l cnp lnp loc 

4 36 0231 2 2 1 i 1 l 1 o 1 0 0 000 

4 39 0234 OOQqOOOOOOO 000 

5 3 0237 0-1 0 o 1 3 0 0 0 0 0 453 

3 6 0242 221illlolOO 000 

5 9 0245 221ilil . 0 100 000 

5 12 0250 221x1110100 000 
5 15 0253 221'illlolOO 000 
|5 16 Q256 221illlo. 100 000 
5 21 026i 2 2 i,i 1 1 i 0 1 0 0 000 
5 24 0264 221ilil 0100 000 
5 27 0267 221illlo 1 0 0 000 
5 30 0272 221ilil0l00 000 
5 33. 0275 221iiliolOQ 000 
5 36 q3oo 221'illlolOO 000 
5 39 q3q3 OOOqqOOOOOO 000 


I 

I 
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PURE ELAS. PROS, USING OLD PCF FQR COMPARISON WITH SOL.N FROM ST£EP32 
CYCLE 2 0 TIME- 0-0000 DT= 1.2700-02 STARTR PROCESSOR . 
STRING 1 POINTS 2 

I J LI LJ X Y > L EP 

5 A 5 3 6.9999993-02 ' 9.9999999-04 1 

3 A 2 3 0,0000000 9,9999999-04 2 

STRING 2 POINTS 2 

I J LI LJ X Y. L EP 

3 3a 2 39 0,0000000 1.1112000+00 3 

3 3« 5 3;J 6.9999996-02 1.1112000+00 4 

FILE 1 DATE 17 MAR 71 TOD 13 Jill 01 BEING WRITTEN ON TAPE 8 
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SHOCK HYDRODYNAMICS, INC, STRIPE CCqe 17 MAR 71 131.11101 


delep 


DELEP 





o o o o o o o o o o o o n o o o o o o n o o o n noon o o o o o o o o o 


PURE EUS, PR08. USING OLD PCF FOR COMPARISON WITH SOLN FROM STEEP32 
CYCLE* 0 T I M F = O.OOOO DT= 1.27QO-02 STARTR PROCESSOR 



I 

X0M 

DXCS 






1.5000-02 

3. 0000-02 




I J 


Y CM 

DYCS 

AREA 

VOL 

U 

0)1 22 

4 

1 . 56/5-02 

3.1750-02 

9 . 2250-04 

9,2250-04 

0,0000 

n 1 2 3 

u; 

4 . 7625-Q2 

3 . 1 7 5 o - 0 2 

9 . 5250-04 

9.5250-04 

0,0000 

0124 

6 

7 . 9 375-02 

3.1750-02 

9 . 525Q-04 

9 . 5250-04 

0.0000 

u 1 ? 5 

7 

1 • 1112-01. 

3 . 1 750-02 

9 . 5 2 5 0 - 0 4 

9 . 525-0-04 

0*0000 

03 26 

a 

1 . 425b" 0 1 

3 . 1 7 5 0 - Q 2 

9 . 5 2 5 n - Q 4 

9 , 525q-q4 

0 • 0000 

G i 2 7 

9 

1 . 7 4 6 2 " 0 1 

3.1750-02 

9 . 5250-04 

9.5250-04 

0.0000 

Q 1 3 f j 

10 

2 . 063/-Q1 

3.i75n-02 

9,525 0 "04 

9 .5230-04 

o, noon 

0131 

11 

2. 30.1.2-01 

3 .1/50-02 

9.5250-G4 

9 . 5 25 [j-U 4 

0-0000 

0132 

12 

2 . oVoZ-Qi 

3 . 1 V 5 q - q 2 

9.525Q-04 

9 . 3250-(j4 

Of 0000 

0133 

13 

3.Q162-Q1 

3.1/50-02 

9,5250-04 

9 . 5250-04 

a- noon 

0134 

14 

3 . 3337-01 

3.1750-02 

9.5250-04 

9 , 5250-04 

0.0000 

0135 

13 

3. 6313“ Qi 

3.1 7 30-02 

9,525g-Q4 

9 , 5250-04 

0,0000 

0136 

1 5 

3.9667-01 

3 . 1.750-02 

9 . 525O-04 

9,5250-04 

n .oono 

013/ 

17 

4 . 2862-01 

3 .1750-Q2 

9 . 5250-04 

9 . 5250-04 

0,0000 


0140 lo f.6037-0i 3.17SQ-02 

0141 l'v 4,9212-Ql 3.175Q-02 

0142 22 2307-01 3.1/5Q-Q2 

0143 21 3.3262-01 3.i75q-q2 

0144 22 2.G7 3/-01 3. 175.0-02 

0i43 23 o , 1 V' 1 2 ~ 0 1 3.175Q-02 

0146 24 ' o . 5(jc / -Qi 3 . 1 7 5 0 - 0 2 

014/ 2 a § • 3262-Gl 3.1./3Q-02 

0l5y 26 7, 143 7-qi 3 • 3.7 5 q “ n 2 

0 1 5 1 2 7 / . 4 6 1 2 “ 0 1 3 . 1 7 5 Q - 0 2 

0152 23 7,7/1’ 7-qi 3.1753-02 

0153 29 q . 0 V 6 2 - 0 1 3.3.750-02 

0154 30 . § . 4137-01 3,l75o^02 

0153 3i §.7312-01 3, 1750-02 

01.56 32 V , 04U7-01 3 • 1 7 5 0 - Q 2 

0 1 5 7 33 y . 3662-01 3 . 1. 7 5 0 - 0 2 

C 1 6 rj 34 v, 6637-01 3*1750-02 

0161 35 i. QUO 1 + 00 3.1 75Q-02 

0162 33 1 .Diiy + 00 3 . 3.750-02 

Q163 37 1. 0636+00 3, 1/50-02 

0164 33 i .0954 + 00, 3.1/50-02 

LI XCS DXCM 

L 3 3.0000-02 3.0000-02 


L 

L I J 

u- 

YCS 

DYC-M 

AREA 

VOL 

TXX 

L 

0121 

i 

-u .qooo 

3.1750-02 

4 .4625-04 

4 .4625-04 

-0,0000 

L 

0122 

4 

5 . 1750-02 

3 . 1750-02 

9.5250-04 

9.525Q-04 

0,0000 

L 

0123 

■ 5 

6 , 3300-02 

3.1/50-02 

9 . 5250-04 

9.5250-04 

Of 0000 

L 

0124 

u 

9 . 523Q-02 

3.1/50-02 

9,5250-04 

9.5250-04 

0. 0000 

L 

0125 

7 

I , 2700-01 

3 . 1 7 5 0 “ 0 2 

9.525Q-Q4 

9,5250-04 

0.0000 

L 

0126 

d 

1 . 5673-01 

3.1/50-02 

9 , 5 2 5 O - 0 4 

9,5250-04 

0,0000 

L 

0127 

9 

1. 9050-01 

3.175Q-02 

9 . 5250-Q4 

9,5250-04 

0,0000 

L 

' 01 30 

10 

2 . 2225*01 

3.1750-02 

9.5250-04 

9 , 5250-04 

0,0000 

L 

0131 

11 

2,3400-01 

3.1/50-O2 ' 

V ,5250-04 

9.5250-04 

0.0000 

L 

0132 

12 

2,8575-01 

3.1750-02 

9 . 5250-04 

9,5250-04 

o.oooo 

L 

0133 

13 

3.175(1-01 

3.1750-02 

9.5250-04 

9,5250-04 

a, oono 

L 

Q 1 3 4 

14 

3 , 4V25-Q1 

3 , 1750-02 

9 , 5250-04 

9,5250-04 

0,0000 

L 

0135 

15 

3 . 8100-01 

3.1/50-02 

9,5250-04 

9 , 3250-04 

o.oooo 

L 

0136 

10 

4.1275*01 

3.17 50-0,2 ■ 

9.5250-04 

9.5250-04 

0,0000 


9 , 5Z50-Q4 

9 ,5250-04 

0,0000 

9,5250-04 

9.5250-04 

0 . OOOQ 

9 , 5230-04 

9 . 5250-04 

0*0000 

9 . 5 2 5 Q - Q 4 

9 ,5250-04 

0,0000 

9 . 5250-04 

9 ,5250-04 

0 . 0000 

9 , 5250-04 

9 . 5230-04 

0.0000 

9 , 5250-04 

9.5250-04 

o , oooo 

9 , 5250-04 

9.5250-04 . 

0-0000 

9 , 525q“ 0 4 

9 • 5 2 5 Q “ 0 4 

0’ 0000 

9,5250-04 

9 >5250-04 

o.oooo 

9 , 525Q-04 

9.525Q-04 

0-0000 

9.5250-Q4 

9.525Q-04 

0,0000 

9 , 525n-04 

9 ,5250-04 

a,oooo 

9.5250-04 

9.5250-04 

0,0000 

9 , 5250-Q4 

9 • 5250*04 

0 * OOOQ 

9 , 5250-04 

9 , 5250-04 

0.0000 

9 . 525Q-04 

9.5250-04 

0,0000 

9 , 5250-04 

9 , 5 2 5 Q 0 4 

0 , 0000 

9.5250-04 

9,5250-04 

0,0000 

9 . 5250-04 

9.5250-04 

0.0000 

9,5100-04 

9 .5100-04 

0,0000 


SHOCK HYDRODYNAMICS, INC, STRIPE COqe 17 MAR 71 13 { 11 * 02 


V 

E 


M 

RHO 

7,5000-01 

0, 

oooo 

2,4926-03 

2 ,7020 + 00 

7 , 5000-01 

0. 

oooo 

2,5737-05 

2,7020+00 

7 .5000-01 

0, 

oooo 

2,5737-03 

2,7020+00 

7 ?5oqq-o1 

o, 

oooo 

2,5737-03 

2. 7q2q+00 

7 • 3000-01 

0 , 

OOOQ 

2,5737-03 

2-7020+00 

7,5000-01 

0, 

oooo 

2,5737-03 

2,7020+00 

7,5000-01 

0. 

0000 

2,5737-05 

2,7020+00 

7 , 5000-01 

0, 

oooo 

2,5737-03 

2. 7Q2Q+00 

7 *5000-01 

0, 

oooo 

2,5737-03 

2.7n2o+00 

7. 5000^01 

0, 

oooo 

2.5737-03 

2,7n2o+00 

0,0000 

0, 

oooo 

2,5737-03 

2* 7060 + 00 

0,0000 

0, 

oooo 

2,5737-03 

2.7Q20+00 

0,0000 

0, 

oooo 

2,5737-03 

2,7020+00 

ai oono 

0, 

oooo 

2,5737-03 

2,7020+00 

0,0000 

0. 

oooo 

2.5737-03 

2.702Q+OO 

n> no-no 

0, 

oooo 

2,5737-03 

2, 7020+00 

n , oooo 

0, 

oooo 

2,5737-03 

2,7020+00 

0 ,0000 

o. 

oono 

2.5737-03 

2.7020+00 

o , n o 0 0 

0, 

oooo 

2,5737-03 

2. 7020 + 00 

n . oooo 

0, 

oooo 

2,5737-03 

2.7020+00 

n • oooo 

0, 

oooo 

2,5737-03 

2*7020+00 

0*0000 

0. 

oooo 

2,5737-03 

2,7n20+00 

0 * 0000 

□ . 

oooo ■ 

2,5737-03 

2*7020+00 

0 ,0000 

0, 

oooo 

2,5737-03 

2 *7020 + 00 

- 0*0000 

D. 

oooo 

2,5737-03 

2, 7 02Q+00 

0,0000 

a, 

ODOO 

2,5737-03 

2 . 7020 + 00 

0,0000 

0 . 

oooo 

2,5737-03 

2.7020+00 

0 ,0000 

0, 

oooo 

2,5737-03 

2.7020+00 

0 • 0000 

-0, 

QOOO 

2.5737-03 

2*7020+00 

0,0000 

0, 

oooo 

2.5737-03 

2,7020+00 

0,0000 

0, 

oooo 

2,5737-03 

2.7020+00 

a , oooo 

0, 

oooo 

2,5737-03 

2. 7Q20 + OO 

■ 0*0000 

o , 

oooo 

2,5737-03 

2 ,7020+00 

a » oooo 

0 , 

oooo 

2,5737-03 

2* 702Q + OO 

0 , 0000 

c, 

oooo 

2,5696-03 

2-7020+00 


TVy 

TZ^ 

TXY 

P + G 

-0 , 0000 

-0,0000 

0,0000 

0.0000 

n. oono 

0,0000 

- 0,0000 

0.0000 

o.oono 

0 . oooo 

0,0000 

o . oooo 

o , nooo 

0,0000 

0,0000 

0.0000 

o , oooo 

o , oooo 

o , oooo 

0.0000 

a. nooo 

o.oooo 

0,0000 

0.0000 

0*0000 

0,0000 

0,0000 

o.ooon 

n s n o o o 

o , oooo 

0 , oooo 

0-0000 

0 , oooo 

0,0000 • 

0,0000 

0 .0000 

n • a ooo 

0,0000 

0,0000 

0,0000 

o , nooo 

0 , QOOO 

0,0000 

0 ,0000 

0,0000 

0 , oooo 

o , oono 

0,0000 

0*0000 

0,0000 

0,0000 

0,0000 

o . n o n o 

0,0000 

0,0000 

0. oooo 


SECTION 7 


SAMPLE PROBLEM: 1-D ALUMINUM ON ALUMINUM IMPACT 



Ax = .03 

The sample problem is a 1-D impact of an aluminum projectile 
l--inch thick on an aluminum target. The yield was set to a very large 
number to simulate a purely elastic material. The equation of state 
TBFT02 was used. 

The data on the next pages are for initial conditions, a 10 cycle 
Integration, and an edit of the velocity vectors at cycle 0. Following the 
data for the edit a graph produced by the SC4020 appears. 


169 




7.1 DATA FOR THE INITIAL CONDITIONS OF THE SAMPLE PROBLEM. 


SIMITAL 

MftXCYC=100» NPAUSE=10» NCOL=3, NROw=52, GEOM=l, 

LBLEFT = 2 , LflRTTE = ?V L9L0tfR=l * LBUPPR=1, 

BTN=8, B0T = 8 * NEWHED= . TRUE . * 

V0WMIS=. TRUE. * SIN=15* S0<J=i6, 

S I N = 1 7 

TMaX=100*, DT = . 0 1 27 » DTMIN=.002, 

S 1 = . 1 6 , Sp=.35, R 1 = 1 • 5 » R?=.5, 

COEVIS=. 3* 

T'/fOMU= . 56 * 

YTELD=1.E+10» 

NOX = 1 , N0Y=2 , X = 0 . , • 09 » OX=.03, 

NR0W = 37 » 

Y=-. 03175*1. 1271, 

NRpG = 2 * XRRE6=.09* .09* YLRFG=0 . * . 3 1 75 * - . 03 1 75 , 

YUPEGr. 31750. 1.5875, 

RHOREG=2*2.702. 

PCON = 2*71 ,4*0. . 1.1867, .7630, 3.445, 1.545, . 9643 ,. 4338 ». 5487 , 1 . 5 , 
PCON(I) =2.702* 
l-it=i* lpoint=i* 

DY=. 03175*. 03175* . . 

VRfG= . 75 , 

YL.REG = 0* , *317, 1.5875, YURFG = . 3 1 8 , 1 . 5875 , 

SEND 

purely elastic test problem for steepbp initial conditions 
startr 

sbdstrg 

NSTRGsl, 

NSTRG=2, 

XSTRG=0. , . 09* 

XSTRG=.09,0. * 

YSTRG=1 .E-20, 1 .E-20, MSTPG = 2. 

YSTRG=2*. 001 ,2*1.1112, 

XSTRG (3) =0. * .09, 

MSTRG ( 2) =2 , 

SEN!) 


4 
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7.2 DATA FOR INTEGRATION OF THE SAMPLE PROBLEM. 


8 9 16 17 l 

JDATAL2 

DT=.005, DTMINs.OO?* NPAUSE=10, MaXCYC=10« TW0MU=.56, 
NEWHED-. FALSE. . COEVIS=0.. 

PCON=2.702»4*o. , 1 • 1867 ,0.763 « 3. *»45, 1.545 *0.9643 » 0.4338, 0.5487 

1 • 5 « 

VONM I S= • TRUE . . SF=. FALSE.. NIM=0, 

SENO 


7.3 data for edit OF SAMPLE PROBLEM. 


SINITAL 

B I M=8 i CYCLE=0» FRESH= .FALSE . » SlN=16» SOU=l7, 

' SEMO • - - • 

EDIT 

JSCA020 

PAPERT=1.* XMIN=0.» XMAX=.OR« YM IN=0 • » YMAX=1.0654, OELX=. 03175, 

DELY=. 03175, FCLASS=2,XSCALE=1 . ,TYPLEN=?3.6,YSCaLE=1 . , 
YMTN=.22225»YMAX=.41275*OELX=.03, 

fmin=o. ,fsc ale= i * » 

LIMEDT = .FALSE. ,plot=. true. ,MGRAPH=1 ,nedit = i ,ooall=. false. , 

P0INTR=1 »NEWHEO=. false. *LAST=. TRUE. , 

REW=.TRUE. ,TYPLEM=.7S»FSCALE = ?5. , 

?ENO 

VELOCITY FIELD (CM/USEC) . 

Z CM RADIUS CM 

$RST ART 
SEND 

CUTOFF 

IffliR 
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CLASTIC TEST PAG6LCM TCR STEEPS* INITIAL CONC I T IONS 


4.IIII10 


' -01 

S.C10X10 
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SECTION 8 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


ARCM (LJ) 
aRCS(J) 
AREAC <L) 
AREAL (L) 
BIN 
BOT 
C ( I J) 
COEVIS 
CMKNEN ( I ) 
CMMGMM ( 1 ) 
CMMSEN ( I ) 
CMSPEN ( I ) 
CMXMOM ( I ) 
CMYMOM ( I ) 
OEM (LJ) 
DEV (LJ) 
DMH (LJ) 
OMV (LJ) 
OUH (LJ) 
OUTOL 

DUV (LJ) 
DVH (LJ) 
DVTOL 

OVV (LJ) 
DXCM (LI ) 
DXCS ( I ) 
DXTOL 

OYCM (LJ) 
OYCS ( J) 

oytol 

E ( I J ) 

EEST ( I J) 
ETOL 

ETOT ( T J) 
FHH ( I ) 

FHV ( J) 

FTT 
FVH ( I ) 

FVV CJ) 

GEOM 

I 

I J 

I MAX 
I M I N 
10 
J 

JMAX 

JMIN 

JO 

JSEN 

KINENG ( I ) 


GLOSSARY OF VARIABLE NAMES 


CUMULATIVE 

CUMULATIVE 

cumulative 

cumulative 

CUMULATIVE 
HOPIZONT *L 


QUADRATIC 


0 


SECTION 

I 

AT 

XEOIT 

SECTION 

I 

AT 

XEDIT 

SECTION 

I 

AT 

XEDIT 

SECTION 

I 

AT 

XEUI T 

SECTION 

I 

AT 

XEDIT 

SECTION 

I 

AT 

XEDIT 


AREA STRESS CELL SIDE LJ 
AREA CELL SIOE J 
AREA BOUNDARY CELL' 

AREA BOUNDARY STRESS CELL 
RESTART INPUT TAPf 
RESTART OUTPUT Tare 
CELL CONTROL WORD 
COEFFICIENT OF VISCOSITY EOP 
CUMULATIVE KTN ENERGY THROUGH 
MOW MAGNITUDE THRU 
MASS * energy THRU 
SPCIFC ENERGY THRU 
X-MOMENTUM THROUGH 
Y-MOMENTUM THROUGH 
SURFACE ENERGY TRANSPORT 
VERTICAL SURFACE ENERGY TRANSPORT 
HORIZONTAL SURFACE MASS TRANSPORT 
VERTICAL SURFACE MASS TRANSPORT 
HORIZONTAL SURFACE U MOMENTUM TRANSPORT 
MINIMUM ABSOLUTE CHANGE IN U VELOCITY 
USED BY PHASEA. 

VERTICAL SURFACE U MOMENTUM TRANSPORT 
HORIZONTAL SURFACE V MOMENTUM TRANSPORT 
MINIMUM ABSOLUTE CHANGE IN V VELOCITY 
USED BY PHASEA. 

VERTICAL SURFACE V MOMENTUM TRANSPORT 
XCM ( 1 ♦ I ) - X C Ml ( T ) 

XCS (LI ) -XCS (LT- 1 ) 

MINIMUM ABSOLUTE POINT DISPLACEMENT IN X DIRECTION 
USED IN RDMOVE. 

YCM ( J* 1 >-YCM( J> 


YCS (LJ)-YCS (LJ-1 ) 

MINIMUM ABSOLUTE POINT DISPLACEMENT IN Y DIRECTION 
USED BY BDMOVF. 

SPECIFIC INTERNAL ENERGY 
ESTIMATED TOTAL ENEPGY/MASS 
MINIMUM ABSOLUTE F 
USED BY MASFLO. 

ESTIMATED TOTAL ENERGY 
HORIZONTAL FORCE ON HORIZONTAL FACE 
HORIZONTAL FORCE ON VERTICAL FACE 
ROTATIONAL FORCE 

VERTICAL FORCF ON HORIZONTAL FACE 
VEPTlCAL FORCE ON VERTICAL FACE 
GEOMETRY. I=SLA8,. 2=CYLIN0RICAL 
cell COLUMN 
cell SUBSCRIPT 
MAXIMUM CELL COLUMN 


MINIMUM CELL COLUMN 
CELL COLUMN OFFSET 
CELL ROW 

MAXIMUM CELL ROW 

MINIMUM CELL ROW 

CELL ROW OFFSET 

SENSE Sw I THE5 0-OFF 1-ON 

KINETIC ENERGY Thru SECTION I AT XE.OIT 
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LBLEFT 

lblowr 

lbkite 

LBUPPR 
LE1T (K) 

LI 

LIJ 

LIMAX 

LIMIN 

L 1 0 

LJ 

LJMAX 
LJMIN 
L JO 

LPOINT <K) 
MASS ( I J) 
MASENG(I) 
MaGMOM ( I ) 
NCOL 
NROW 
PCON (L) 

POS 

PTOL 

PA (L»M*N) 
RB <L*M, N) 
RHOTOL 

O 

SMASS(I) 
STT A (LIJ) 
STTEST (LIJ 
STTH 
STTV 

S X X A (LIJ) 
SXXEST (LIJ 
SXXH 
SXXV 

S Y V A (LIJ) 
SYYEST (LIJ 
SYYH 
SYYV 

TTT (LIJ) 
TTTTOL 
TXX (LIJ) 
TXXTOL 
T X Y (LIJ) 

T X Y A (LIJ) 
TXYEST (LIJ 
TXYH 
TXYTOL 
TXYV 

TYY (LIJ) 
TYYTQL 
TWOMU 
U(IJ) 

UAVE ( I J) 
UEST ( I J) 
USKING ( I ) 
USMASS ( I ) 
USMGMM ( I ) 
USMSEN ( 1 ) 


BOUNDARY i=free*2=fixed 
BOUNDARY i=free*2=fixeq 
BOUNDARY 1=FREE»2=FIXE0 
BOUNDARY 1 =FREE « 2=F I XEO 
OF FIT FOR EQUATION OF STATE 


COLUMN 

COLUMN 

offset 


LEFT 
LOV/R 
RITE 
UPPR 
TYPE 

LATTICE COLUMN 
LATTICE 
MAXIMUM 
MINIMUM 
LATTICE 
LATTICE 
MAXIMUM 
MINIMUM 
LATTICE 
POINTER 
CELL MASS 

MASS*ENERGY THROUGH 
MOMENTUM mag thru 
number OF CELL COLUMNS 
NUMBER OF CELL ROWS 
EQUATION OF STATE CONSTANTS 


LATTICE 
LATTICE 
COLUMN 
ROW 

LATTICE ROW 
LATTICE ROW 
ROW OFFSET 
TO LOCAT ION 


OF EQUATION OF STATE CONSTANTS 


SECTION 

SECTION 


AT 

AT 


xeqit 

XEDIT 


TRUE IF 
MINIMUM 
MINIMUM 
MAXIMUM 
MINIMUM 
USEO BY 
MASS 
S/TT 
ESTIMATED 


ANY REAL STRING POINTS 
ABSOLUTE PRESSURE 
cut point of segment 
cut point of segment 
ALGEBRAIC DENSITY (RHO=0. 
80*0 VF. 

through section i at xedit 
* A ( N ) 


in problem 


IF LESS THAN RHOTOL) 


s/tt 


STRESS 

flow 

STT 

ACROSS 

HORIZONTAL 

SURFACE 

STRESS 

FLOW 

STT 

ACROSS 

VERTICAL 

surface 

S/X.X * A (N) 







ESTIMATED S/ Y X 






stress 

FLOW 

sxx 

ACROSS 

HORIZON! 

AL 

SURFACE 

STRESS 

FLOW 

sxx 

ACROSS 

VERTICAL 

SURFACE 

S/YY * A (N) 







ESTIMATED S/YY 






stress 

FLOW 

SYY 

ACROSS 

HORIZONT 

AL 

SURFACE 

stress 

FLOW 

SXY 

ACROSS 

VERTICAL 

SURFACE 


SIGMA/THETa-THETa 
MINIMUM ABSOLUTE TTT 
SIGMA/XX 

MINIMUM ABSOLUTE TXX 

SIGMA/XY 

T/XY * A ( N ) 

ESTIMATED T/XY 

STRESS FLOW TXY 
MINIMUM ABSOLUTE TXY 
STRESS FLOW TXY 

sioma/yy 

MINIMUM ABSOLUTE TYY 

?*MU for material 
X VELOCITY 
X VELOCITY AVERAGE 


ACROSS HORIZONTAL SURFACE 
ACROSS VERTICAL surface 


X VELOC 

TTY E 

ST I MATE 



KIN ENG 

PER 

steradian 

AT 

XEDIT 

MASS 

PER 

STERAD I AN 

AT 

XEDIT 

MOM MAG 

PER 

steradian 

AT 

xedit 

mas*eng 

PER 

steradian 

AT 

XEDJT 


USXMOM ( I ) 
USYMOM ( I ) 
UTOL 

V(IJ) 
VAVE(IJ) 
VEST ( I J) 
VOLC (L) 
VOLTOL 
VOLL (L) 
VTOL 

XCM ( I ) 

XCS (LJ) 
XMOM(I) 
YMOM ( I ) 
YCM ( J) 

YCS (LJ) 
YIELD 


X-MOM PER STERADIAN AT XEDIT 

Y-MOM PER STERADIAN AT XEDIT 

MINIMUM ABSOLUTE U 
USED IN MASFLO. 

Y VELOCITY 

Y VELOCITY AVERAGE 

Y VELOCITY ESTIMATE 
VOLUME BOUNDARY CELL 

MINIMUM VOLUME RATIO FOR ACCELERATION 
VOLUME BOUNDARY STRESS CELL 
MINIMUM ABSOLUTE V 
USED IN MASFLO AND PHASEA. 

X center OF MASS 
X CELL SIDE 

X-MOMENTUM THROUGH SECTION I AT XEDIT 
Y-MOMENTUM through SECTION I AT XEOir 

Y CENTER OF MASS 

Y CELL SIDE 
YIELD CONSTANT 


SECTION 9 


SUMMARY AND RECOMMENDATIONS 


Under the current contract, the STEEP code has been converted to 
the Marshall Space Flight Center EXEC VIII system. The converted code 
has been designated STEEP32. This document describes STEEP32 from the 
User's point of view. Each major step in a STEEP32 solution is illustrated 
in a sample problem. There is a detailed discussion of the internal 
organization of the code. 

To facilitate the utilization of STEEP32 by MSFC personnel it is 
recommended that workshop sessions be conducted. The utilization of 
the STEEP32 documentation to set up and run problems of current interest 
would be included in these sessions. 

It is further recommended that provisions be made for continuing 
technical staff services. Such services would include additional 
instruction and could be extended to programming assistance in making 
code modifications. 

Finally, additional consideration should be given to improved 
numerical techniques which are now available. In particular, the basic 
formulation of STEEP32 has been extended to multiple materials. This code 
is called SHESAM (an acronym for Shock Hydrodynamics Eulerian with 
.Strength and Multiple materials) . Reference 5 describes the results of an 
application of SHESAM to the study of debris formed by a shape charge 
perforation . 
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APPENDIX A 


EQUATIONS OF STATE 

Included in the following are two of the most frequently used equations 
of state. Additional equations of state which have been used in conjunction 
with codes developed by Shock Hydrodynamics Incorporated can be found in 
Reference 4. 

EQUATION OF STATE SUBPROGRAMS SPECIFICATIONS 
Nomenclature Units* 


p 

material density 

gm/cc 

e 

specific internal energy 

, rt l2 z Mb-cc 

10 ergs/gm, 

^ ^ gm 

P 

pressure 

megabars 

c 

adiabatic sound speed 

cm/p s 

o 


♦with certain exceptions 

Definitions 



p o 

normal material density 


e 

o 

rest energy 


e' = e + e 

o 



T1 = P/P Q 

relative density 


> =11-1 

compression 


V = 1/p 

specific volume 


e = p e 
o 

energy density 



Universal Constants (C (I) is first for particular fit) 

C(I) p Q normal material density 

C(I + 1) c min a) the minimum sound speed 


*C(I + 2) | p| min If (| p | < C(I + 2))p is replaced by 0; 

T) is replaced by 1 
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Universal Constants (continued) 


*C(I + 3) 

e min 

*C(I + 4) 

P min 

Inputs 

T|, e 

Outputs 

P, c^ where 


If (e < C(I + 3)), e is replaced by 0 
If (P < C(I + 4)), P is replaced by 0 
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FIT 1: LASL Metal Fit 


FORM: 


P — A + Be 4- Ce' 

e + e 
o 


A = n (a x + a 2 i n| ) 

B = b Q + b l p + b 2 p 2 
c = c 0 + c^ 


° 2 - p (c * + .) U- 

+ B'e + C e 2 + 

-j (B + 2Ce 

-P)} 




■n 



A' a aj^ + 2a 2 | \x j 





B* = e b 1 + 2b^ 





r—H 

o 

II 

b 




' 

CONSTANTS (13) 

Maqnesium 1 

Aluminum 2 

Iron 3 

Beryllium' 

cm Pq 

1.735 

2.702 

7.86 

1.845 

C (l + 5) a 

0.56645 

1.1867 

7.78 

0.9512 

C(l + 6) a 2 

0.33433 

0.7630 

31.18 

0.3453 

C(l + 7) b Q 

2.2178 

3.455 

9.591 

0.9269 

C (I + 8) b x 

0.87104 

1.545 

15.676 

2.948 

C (I + 9) b 2 

0.48136 

0.9643 

4.634 

0.5080 

C(I+ 10) c Q 

0.41626 

0.4338 

0.3984 

0.5644 

C(I+11) c 1 

0.58904 

0.5487 

0.5306 

0.6204 

C (I + 12) e Q 

1.5 

1.5 

9.0 

0.8 


1 Source for MG LASL (R. K. Osborne, Letter Dated 12 Oct. 1964). 

2 Source of data: L-401, 8 Jan. 1960, R. L. Bjork, The RAND Corporation . 
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FIT 12: Tabular Fit 2 (Shock Hydrodynamics) 

FORM : Double interpolation /extrapolation on log r\ , log e' (e' = e + e ) 

in table of P(log T] , log e') ° 



2 


_L 3P + 3P 3_e 




p Q 3 t) 3 e 3 r\ 


3P 


P(T1 i+ 1' e ’ } - p ^i' 


ST] 


^i+l^i 


IP 


P(ri, e! + x ) - P (ri, e!) 


3 e 



when 

i 

is 

such that T] £ r} s r). ^ 


j 

is 

such that ej £ e ^ e j +1 


and extrapolation is done in the usual fashion. 

CONSTANTS 

(12 + NETA + NETA * NE) 

C(I) 

p o 

C(I +5) 

a If a = 0 and P < P min P is replaced by 0 
If a > 0 and P < P . , P is replaced by P . 

C (I + 6) 

date, date of table generation 

C(I +7) 

e Q , free energy 

C (I +8) 

NETA (floating point) number of (log T|)'s 

C (I + 9) 

NE (flating point) number of (log e')'s 

C(I+10) 

min log e' 

C(I +11) 

Alog e' (table is uniform in Alog d') 

C(I + 12) 

log Tj^ (table allows for non-uniform log T]) 

C(I+ 13) 

log r\ 2 
etc. 
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20 NSX AJVjs 

C (I+NETA + u ) 
0(1 +XETA + i 2 ) 


(12 * neta + nexa 


109 \ 


ax 


P (log 


NE) 


C(J+ «m + 13) V 

(log V log 

Cl, etc - 

0(1 +W£TA + n 

1 +N ™*m )P(l0gri 

neta • J °g 


o 


(contj nuec jj 
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APPENDIX B 


FORMAT OF STEEP 3 2 RESTART TAPES 


A STEEP32 restart tape consists of a sequence of logical files. 
Each logical file corresponds to a particular computational cycle that 
the STEEP user has elected to save. 

Each logical file consists in turn of a sequence of logical records. 
Each logical record is defined by a record type. Record type 1 is a 
header record containing the file number of the file, the time and date 
it was created and the record type. All succeeding record types begin 
with the same format as record type 1. 

Record type 2 contains the variables from common blocks CONSTS 
and MATBLK. Record type 3 contains common block BOUNDY. There are 
multiple type 4 records. The first record type 4 contains the cell 
dictionary. Each succeeding record type 4 contains the cell variables 
for a particular column. These correspond to the first twelve elements 
of the drum file. Record type 5 contains CUTBLK as well as the total 
number of allowable elements in any one of the arrays in CUTBLK. 

Record type 6 contains the grid coordinates XCS and YCS. There 
currently is no record type 7. Record type 8 is a trailer record for the 
file. Record type 9 is a logical end of tape record. 

' The STEEP32 restart tape is written using FORTRAN binary write 
statement. 


o 
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